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A  method  For  the  construction  of  a  rectangular 
geometric  dual  from  a  Deltahedron  based  maximally  planar 
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CHAPTER  1 


INTRODUCTION 

The  problem  of  whBre  to  locate  different 
facilities  uiithin  a  structure  is  a  very  old  one  indeed. 
Whenever  a  building  serves  more  than  one  function  with 
each  function  having  specific  equipment  or  space 
requirements,  choices  must  be  made  to  determine  the  best 
location  for  each  function.  Even  thB  simple  problem  of 
locating  a  bed,  fireplace,  and  table  within  a  cabin 
requires  choice  among  differing  alternatives.  This 
problem  however,  is  not  limited  tD  location  of  rooms  or 
functions  within  a  building.  Extensions  can  be  made  to 
include  problems  ranging  from  the  location  of  different 
buildings  on  a  site  to  electronic  components  on  a  circuit 
board.  Many  approaches  to  this  problem  have  been  taken 
over  a  great  span  of  time.  One  approach  sometimes 
referred  to  as  iconic,  includes  building  models  of  the 
different  components  and  physically  placing  them  in 
different  locations  within  a  model  of  the  building.  The 
analog  approach  is  one  that  transforms  the  original 
problem  into  some  analogous  problem  and  then  solves  this 
analog  problem.  The  solution  for  the  original  problem  is 
then  obtained  by  a  reverse  transformation.  The  approach 


that  as  of  late  has  had  by  far  the  most  attention  is  the 
symbolic  or  mathematical  approach. 

This  thesis  deals  with  the  extension  of  several 
specific  mathematical  approaches.  In  particular,  the 
development  of  the  spacial  relationships  infered  by  the 
results  of  a  special  class  of  graph  theoretic  methods 
known  as  Oeltahedron  Heuristics. 

The  purpose  of  this  thesis  is  to  develop  a 
systematic  approach  to  construct  a  rectangular  geometric 
dual  from  these  Deltahedron  based  adjacency  graphs  and 
include  areas  to  form  a  block  plan.  Chapter  2  describes 
the  problem  as  well  as  some  past  work  in  the  area.  In 
addition  to  a  systematic  approach  for  developing  a 
rectangular  geometric  dual  and  its  block  plan,  a  computer 
implementation  of  this  method  is  included  in  chapter  3. 
Comparisons  with  two  other  computerized  methods  are  given 
in  chapter  4  while  chapter  5  contains  conclusions  and 
suggestions  for  further  work. 


CHAPTER  2 


PROBLEM  STATEMENT  AND  PAST  WORK  IN  THE  AREA 

The  general  purpose  of  all  of  the  layout  methods 
proposed  Is  to  specify  locational  relationships  between 
Facilities  so  as  to  optimize  some  performance  criterion. 
These  relationships  are  generally  of  two  forms,  the 
adjacency  of  facilities  and  the  distance  between 
facilities.  The  most  common  objective  functions  used  to 
measure  the  performance  criterion  are  maximization  of 
total  achievable  adjacencies  and  minimization  of  total 
transportation  cost.  When  maximizing  the  sum  of 
adjacencies,  each  adjacency  between  two  facilities  has 
some  specified  score  and  the  total  of  all  adjacencies 
realized  represents  this  total  adjacency  score.  The 
minimization  of  total  transportation  cost  usually  assumes 
that  transportation  cost  is  a  function  of  distance  and 
therefore  the  overall  pairwise  distance  between 
facilities  that  have  some  material  being  transferred  must 
be  minimized. 

S=l_Claaaical-LawQut-6BBCQQCbea 
The  first  formal  mathematical  model  of  the 
facility  layout  problem  was  in  the  form  of  the  Quadratic 
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Assignment  Problem  proposed  by  Koopmons  and  Beckmann 
C 19573.  This  formulation  takes  the  approach  of  dividing 
each  facility  into  somB  number  of  equal  size 
subfacilities,  usually  using  the  size  of  the  smallest 
facility.  The  task  is  then  to  assign  each  subfacility  to 
a  location  on  an  orthagonal  grid  representing  the  planar 
site,  so  that  the  total  transportation  cost  is  minimized 
and  that  each  facility  occupies  a  contiguous  region.  It 
has  been  shoun  that  this  problem  has  no  algorithm  for  its 
solution  that  is  polynomial ly  bounded  in  time  and  belongs 
to  the  class  of  problems  termed  NP  complete.  This  means 
that  only  relatively  small  problems  can  be  solved  to 
optimality  using  this  method.  Therefore,  attempts  have 
been  made  to  find  a  good  heuristic  to  provide  solutions 
to  this  problem.  Some  of  the  well  known  methods  are 
briefly  described  below. 

2-1.1  Terminology,  Notation,  and  Definitions 

The  following  terms  and  notation  are  defined  in 
the  context  of  facility  layout. 

Lll-CQDatCUCtian.UBUCiatiC.  a  construction  type 
heuristic  is  one  that  constructs  a  layout  by  adding 
facilities  one  at  a  time  until  a  completed  layout  is 
achieved . 

LSl.lODCauoOBDt-UeuciatlC*  An  improvement 
heuristic  is  one  that  requires  an  initial  layout  as 


input.  The  heuristic  then  improves  the  layout  by  some 
local  exchange  technique  until  no  further  improvements 
can  be  made . 

L31_KBlatiQDSbiB_Cbact •  The  relationship  chart, 
or  REL  chart,  is  an  attempt  to  quantify  the  importance 
of  relationships  between  facilities  using  closness 
ratings  Ctluther,  19611 .  The  closeness  rating  is  a  score, 
Rjj ,  that  is  achieved  when  the  two  appropriate 
facilities  are  adjacent.  The  ratings,  their  definitions, 
and  frequently  used  scores  for  two  common  methods  are 
listed  in  Table  2.1. 

Lil-fidiacency •  Generally  two  facilities  are 
considered  adjacent  if  they  share  a  common  wall  or 
divider  of  some  minimal  tolerance  length  that  separates 
one  from  the  other.  One  exception  to  this  definition  is 
the  criterion  of  ALDEP  which  in  addition  to  the  above 
description,  considers  two  facilities  adjacent  if  they 
are  diagonal  to  one  another  at  the  meeting  of  four  walls. 

LSl-lbitial-Layaut •  The  initial  layout  is  the 
layout  used  for  a  starting  point  in  improvement  type 
heuristics . 

£.6  2  _E  law -Dot  a  •  This  is  a  matrix,  sometimes 
refered  to  as  a  From-To  chart,  that  represents  the 
number  of  trips  or  volume  of  material  flow  per  time 
period  being  made  from  one  facility  to  another. 


I2J_Coat_DatO .  This  is  also  a  matrix  however  it 
contains  the  cost  to  move  one  unit  of  distance  between 
each  facility. 

LSl-Elaat-Lauaut •  Since  the  majority  of  layout 
planning  has  dealt  with  the  design  of  manufacturing 
structures,  the  building  or  collection  of  buildings  is 
commonly  referred  to  as  the  plant;  hence  the  term  plant 
layout . 


Table  5.1  Common  REL  Chart  Ratings,  Definitions, 
and  Scores 


Rating 


Definition 


A 

E 

I 

□ 

U 

X 


Absolutely  necessary 
Especially  important 
Important 

Ordinary  closeness  OK 

Unimportant 

Undesirable 


Score 

._dLDEE _ CQBELQE. 

64  6 

16  5 

4  4 

1  3 

0  5 

-1054  1 


5-1.5  Muther’s  Systematic  Layout  Planning 

Muther,  C 19611  developed  the  organized  approach 
to  plant  layout  known  as  Systematic  Layout  Planning 
C5LPJ.  The  three  main  areas  of  concern  for  this  method 
are  Analysis,  Search,  and  Selection  as  illustrated  in  the 
method  schematic  shown  in  figure  5.1. 


-<SV: 
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within  the  plant  is  collected  in  the  form  of  a  flew  and  a 
cost  from-to  chart.  Additional ly ,  quantifiable 
information  about  the  desirability  of  having  each  pair  of 
facilities  within  the  plant  adjacent  to  one  another  is 
collected  in  the  form  of  a  REL  chart  Csee  figure  E.SCaJl. 
The  information  from  these  three  is  then  used  to  come  up 
with  a  relationship  diagram.  The  relationship  diagram  is 
constructed  by  arranging  equal  area  squares  that 
represent  each  facility  into  different  configurations 
until  one  is  found  that  has  the  desired  level  of 
preferred  properties  measured  by  the  from-to  and  REL 
charts  Csee  figure  2.2Cb33.  This  is  often  an  iterative 
trial  and  error  scheme  that  ie  performed  manually  with 
evaluation  often  done  very  subjectively  and  therefore 
many  and  possibly  preferred  arrangements  may  be 
overlooked.  Space  requirements  for  each  facility  are 
then  determined  as  well  as  the  total  available  space. 

LSl-Seaccb •  Tbs  search  operation  is  started  by 
developing  several  space  relationship  diagrams  Csee 
figure  2.3Ca3.  This  involves  combining  the  relationship 
diagram  with  the  space  requirements  and  space 
availability  to  construct  diagrams  that  have  the 
relationships  and  areas  suggested  during  the  analysis 
stage.  These  space  relationship  diagrams  are  then 
condensed  into  a  block  plan  as  illustrated  in  figure 


Figure  2.2 — Continued*  Relationship  Diaoram 

2.3Cbl.  This  block  plan  is  finally  combined  with  any 
modifying  considerations  and  practical  limitations  that 
are  developed,  to  come  up  uith  alternatives  for  the  plant 
layout . 

L3 2 -Select ion ■  The  final  operation  is  to  decide 
among  the  alternatives  or  to  make  any  data  changes  that 
prove  necessary  and  repeat  the  process. 

All  other  methods  presented  here  fit  uithin  the 
general  context  of  this  procedure.  Any  layout  will 
involve  collecting  data  and  some  selection  among 
alternatives.  The  difference  arises  uith  the  choice  of 
the  method  one  uses  to  construct  the  block  plan  from  the 


Figure  2.3.  CqD  Space  Relationship  Diagram  CbD  Block 
Plan 


2-1.3  ALDEP 

A  method  that  mas  developed  within  IBM  and 
original lg  presented  bg  SeehoF  and  Evans  C1S673  is  called 
the  Automated  Logout  Design  Program,  commonlg  referred  to 
as  ALDEP.  ALDEP  is  a  construction  tgpe  heuristic  as  it 
requires  no  initial  solution  to  begin,  however  it  uses 
its  past  solutions  as  a  basis  for  comparing  new  ones  to 


IE 


see  if  any  improvement  has  been  made  and  therefore  some 
improvement  does  take  place.  ALDEP  divides  each  facility 
into  subfacilities  of  same  common  square  dimension  based 
upon  thB  scale  specified.  A  facility  is  then  chosen  at 
random  and  layout  is  begun  from  the  upper  left  corner  of 
the  layout.  The  subfacilities  of  the  initial  facility 
are  added  to  the  layout  in  vertical  strips  of  a  specified 
'sweep  width’  until  its  area  is  exhausted.  The  REL  chart 
is  then  scanned  for  a  facility  that  has  an  A  or  E  rating 
with  the  existing  facility  and  it  is  then  placed  in  the 
layout.  As  before  the  new  facility  is  laid  in  a  strip 
fashion  until  its  area  is  exhausted.  The  vertical 
scanning  nature  of  these  strips  is  illustrated  in  figure 
2.4.  This  addition  process  is  then  repeated  until  no 
facilities  remain  or  until  there  are  no  facilities  with 
an  A  or  E  rating  with  the  last  facility  added.  If  the 
latter  is  the  case,  a  facility  is  chosen  at  random  and 
the  process  is  continued.  The  score  for  this  method  is 
found  using  the  values  from  REL  chart.  The  eight  squares 
that  surround  each  facility  are  scanned  and  the  score 
recorded.  After  a  score  is  recorded  it  is  deleted  from 
the  matrix  to  eliminate  the  possibility  of  including  the 
same  adjacency  twice.  The  total  of  these  values  is  the 
score  for  the  layout.  Usually  the  entire  process  is  run 
many  times  with  each  improvement  in  score  becoming  the 
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new  goal  for  the  program  to  attain.  Runs  that  do  not 
achieve  the  goal  are  rejected  and  the  entire  process 
stops  when  no  improvement  is  made.  Alternatively,  a 
collection  of  good  solutions  can  be  developed  to  provide 
different  options  for  the  selection  process.  An  example 
of  the  output  produced  is  included  in  chapter  4. 


Figure  2.4.  (Jertical  scanning  pattern  used  by  ALDEP 

2-1.4  CORELAP 

CORELAP  is  the  acronym  for  Computerized 
Relationship  Layout  Planning  and  uias  developed  by  Lee  and 
Moore  C 1967 3.  A  number  of  improvements  to  the  original 
method  have  been  added  since  its  introduction  and  the 
version  known  as  CORELAP  B  will  be  discussed  here.  As 
with  ALDEP,  this  is  a  construction  type  heuristic.  This 
method  begins  by  choosing  the  first  facility  according  to 
its  Total  Closeness  Rating  CTCR3,  calculated  far  facility 
i  by  summing  the  REL  chart  scores  from  facility  i  to  all 
others.  The  facility  with  the  highest  TCR  is  chosen  to 


be  added  first,  and  placed  in  the  center  of  the  layout. 
Next  a  facility  that  has  an  A  adjacency  score  with  the 
first  facility  is  selected.  If  no  facility  with  an  A 
rating  is  found,  an  E  rating  is  searched  for.  If  no  E 
rating  is  found,  the  method  continues  doun  the  hierarchy 
of  scores  until  a  U  is  reached.  If  no  facility  uith  a 
score  of  U  or  better  is  found,  the  facility  uith  the 
highest  TCR  is  chosen.  If  there  is  more  than  one 
facility  uith  the  same  scare,  the  facility  uith  the 
highest  TCR  is  chosen.  The  same  type  of  search  is 
employed  at  all  sucessive  steps  uith  the  search  started 
by  looking  far  a  facility  with  an  A  adjacency  to  the 
first  facility.  If  none  is  found,  an  A  adjacency  uith 
the  second  facility  is  desired,  folloued  by  an  E  uith  the 
first,  an  E  uith  the  second,  an  I  uith  the  first,  etc. 

All  facilities  are  added  to  the  exterior  of  the  existing 
arrangement  and  are  rectangular  in  shape.  They  are 
placed  in  a  position  that  uill  yield  the  highest 
placement  rating  and  boundary  length, where  the  boundary 
length  is  the  length  of  the  boundaries  common  to  the  new 
facility  and  the  existing  layout.  Some  different 
configurations  possible  are  illustrated  in  figure  2.5. 

The  placement  rating  is  the  sum  of  the  weighted  ratings 
between  the  department  being  added  to  the  layout  and  its 
neighbors  if  it  is  placed  there.  The  weights  are 


assigned  to  the  adjacency  ratings  by  the  user.  Therefore 
the  scorB  used  far  the  TCR  is  not  necessarily  the  same  as 
that  used  to  score  the  placement  of  each  facility  within 
the  layout.  An  example  of  the  output  from  this  method  is 
also  included  in  chapter  4. 


I-* —  3  - 4^2*! 


Figure  2.5.  CORELAP’s  placement  method 

2-1.5  CRAFT 

CRAFT  is  an  improvement  type  heuristic  and  uias 
introduced  by  Armour  and  Buffa  C1963],  In  addition  to 
differing  from  ALDEP  and  CORELAP  in  the  type  of  heuristic 
used  [construction  versus  improvement],  CRAFT  employs  a 
entirely  different  method  for  evaluating  a  layout. 
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Unlike  ALDEP  and  CORELAP,  CRAFT  attempts  to  minimize 
transportation  cost  where  this  cost  is  expressed  in  terms 
of  distance  traveled.  This  is  therefore  an  attempt  to 
provide  a  solution  to  the  DAP  mentioned  earlier.  As  an 
improvement  heuristic,  CRAFT  requires  an  initial  layout 
in  order  to  apply  its  improvements.  The  score  for  a 
layout  is  the  cost  per  unit  distance  Ccost  data!  to  move 
an  item,  multiplied  by  the  rectilinear  distance  between 
facility  centroids,  multiplied  by  the  number  of  trips 
required  Cflow  data!,  for  all  pairs  of  facilities  in  the 
layout.  The  next  step  is  to  consider  the  exchange  of  two 
or  three  facilities  within  the  layout.  The  possible 
combinations  include  13  two-way  interchanges,  S3 
three-way  interchanges,  33  two-way  followed  by  three-way 
interchanges,  43  three-way  followed  by  two-way 
interchanges,  and  53  the  best  of  two-way  and  three-way 
interchanges.  Exchanges  of  facilities  are  only  possible 
if  the  facilities  are  adjacent  to  one  another  or  if  their 
areas  are  equal.  The  search  for  the  best  of  these  is 
done  by  interchanging  the  centroids  which  are  used  in  the 
distance  calculations  as  an  estimate  of  the  actual  cost. 
The  best  exchange,  lowest  score,  is  then  mode  and 
centroids  recalculated  according  to  the  new  shape  of  the 
facilities.  If  a  savings  still  exists  the  process 
continues  and  if  not  the  old  layout  is  maintained  and  a 
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different  interchange  is  attempted.  When  no  improvements 
can  be  made  the  process  stops .  A  drawback  with  the 
method  is  that  there  appears  to  be  no  consistent  method 
for  the  physical  exchange  of  adjacent  facilities  of 
varying  areas. 

SrS-GcQBb-IbsQceticaJL-acQCQQcbea 
2-2.1  Terminology,  Notation,  and  Definitions 

The  following  terminology  and  notation  is 

defined: 


i 

i 


£12  JaCUBb •  A  graph  is  a  pair  of  sets  CO, El  where 
U  is  finite  and  not  empty.  The  elements  of  U  are  called 
vertices  and  the  elements  of  E  are  distinct  pairs  of 
vertices  called  edges.  If  there  is  no  direction 
associated  with  the  edges,  they  are  known  as  undirected 
edges.  If  all  edges  are  undirected,  the  graph  is  said  to 
be  an  undirected  graph. 

£.22 JlleiBbted-ScaBb .  a  graph  that  has  a  weight, 

W_,  assigned  to  each  edge,  e,  is  known  as  a  weighted 
graph  with  W_  usually  being  an  element  of  the  positive 
real  numbers . 

L32_CQ(BBlote_GcaBb •  A  complete  graph,  denoted  K^, 
is  one  in  which  all  pairs  of  vertices  are  joined  by  an 
edge.  A  complete  undirected  or  symmetric  graph  has 
CnCn-1 11/2  edges. 


ISJ-ElDDOC-BrODb •  A  graph  is  said  to  be  planar  IF 
it  can  be  drawn  in  the  plane  such  that  no  two  edges 
intersect  except  at  a  vertex  to  which  both  are  incident. 

I52-D02iDQllU-ElBDBE-GCDQb .  A  graph  is  said  to  be 
maximal ly  planar  IF  it  not  possible  to  add  an  edge  and 
still  maintain  planarity.  Due  to  the  Fact  that  all  Faces 
oF  a  maximally  planar  graph  are  triangles,  a  maximally 
planar  graph  is  oFten  known  as  a  triangulation.  A 
Maximally  planar  graph  contains  3n-6  edges  CEuler,  17523. 

t,61_lBtCBhOClCQn  •  A  tetrahedron  CK^D  is  a  complete 
graph  on  Four  vertices  which  is  also  maximally  planar 
Csee  Figure  2.63. 


Figure  2.6.  Tetrahedron 

£21_0BltabBSlCQn  •  a  deltahedron  is  a  graph  that  is 
constructed  by  beginning  with  a  tetrahedron  and  adding 
vertices  by  the  insertion  oF  an  additional  vertex  into  a 
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triangle  and  adding  edges  From  the  new  vertex  to  each  of 
the  three  vetrices  that  define  the  triangle.  Due  to  this 
fact  a  deltahedron  must  contain  at  least  one  vertex  oF 
degree  three  Cthree  edges  incident  with  itl. 

IB3_0QBimollu_BlDDor_edjQceDCu-GroDb.  a  maximal ly 
planar  adjacency  graph  is  a  maximally  planar  graph  whose 
edges  represent  adjacency  between  pairs  of  facilities. 

ISJ-GBOOBtriC-Oual .  The  geometric  dual  of  the 
maximally  planar  adjacency  graph  is  a  spacial 
representation  of  the  facilities  that  orB  represented  by 
the  vertices  of  the  graph.  The  edges  of  the  graph 
represent  the  adjacency  of  two  facilities  in  the  dual. 

If  a  graph  is  maximally  planar  then  its  dual  is  also 
maximally  planar  or  in  other  words  no  further  adjacencies 
in  the  dual  can  be  established  without  violating  the 
planarity  of  the  dual  C Whitney,  1931 D. 

IlDJ-BBCtQDgulDr_BBQmBtric_DuBl .  For  this 
discussion,  a  rectangular  geometric  dual  is  a  geometric 
dual  that  contains  only  rectangular,  L  and  T  shaped 
areas . 

All  graph  theoretical  approaches  presented  here 
are  of  the  construction  type.  One  starts  with  a  complete 
graph  on  N  vertices  corresponding  to  a  REL  chart  with 
zero  weight  edges  added  if  necessary,  and  attempts  to 
find  a  maximally  planar  subgraph  on  the  complete  graph 
that  has  maximum  weight  since  without  lass  of  generality, 


with  nonnegative  weights,  an  optimal  solution  will  be 
maximally  planar.  The  problem  of  starting  with  the 
complete  graph  and  deleting  edges  until  it  is  maximally 
planar  is  a  relatively  difficult  and  very  time  consuming 
problem  due  to  the  methods  required  to  check  for  maximal 
planarity .  The  methods  shown  here  use  construction 
techniques  that  start  with  either  a  graph  that  is  not 
maximally  planar  and  iteratively  build  it  up  until  it  is 
maximally  planar  or  a  graph  that  is  maximally  planar  and 
then  add  vertices  and  edges  to  it  in  a  specific  manner  so 
that  it  will  always  remain  maximally  planar.  Several  of 
the  methods  start  with  a  complete  graph  on  four  vertices, 
Ki£ .  There  are  basically  two  methods  for  determining 
which  four  vertices  should  make  up  this  initial 
tetrahedron.  The  first  is  the  greedy  approach  which 
finds  the  highest  weight  tetrahedron  among  all 
possibilities.  The  other  is  formed  by  first  summing  the 
scores  of  all  columns  from  the  square  adjacency  matrix. 
The  vertices  are  then  sorted  in  non-increasing  order 
according  to  these  column  sums.  Then  the  vertex  with  the 
highest  adjacency  rating  to  all  other  vertices  is  chosen 
first.  It  has  been  shown  CSiffin,  19B4D  that  there  is 
empirically  no  clear  difference  in  final  triangulation 
solution  quality  for  either  starting  procedure.  The 
objective  of  all  methods  that  follow  Cwith  the  exception 


□f  Super  Deltahedronl  is  to  maximize  the  adjacency  score 
inhere  the  values  of  having  two  facilities  adjacent  are 
the  same  as  those  used  in  ALDEP. 

2-2.2  The  Ulheel  Expansion  Heuristic 

The  Wheel  Expansion  Heuristic  CEades,  Foulds,  and 
Gif fin,  19B21  begins  with  an  initial  tetrahedron  and  uses 
an  operation  known  as  a  wheel  expansion  to  add  sucessive 
vertices  to  the  graph.  It  has  been  shown  that  the  wheel 
expansion  operation  is  sufficient  to  form  all  maximally 
planar  graphs  possible  if  the  starting  point  is 

.  An  example  of  wheel  expansion  is  illustrated  in 
figure  2.7.  The  choice  of  vertex  and  location  for 
expansion  involves  finding  the  vertex  and  expansion  point 
that  has  the  highest  increase  in  adjacency  score. 


Figure  2.7.  Wheel  Expansion 

2-2.3  The  Greedy  Heuristic 

The  idea  behind  the  Greedy  Heuristic  CFoulds, 
Gibbons,  &  Giffin,  19851  is  very  straight  forward. 

First,  all  edges  are  listed  so  that  all  edges  with  A 
values  ars  first  followed  by  those  that  have  a  value  of  E 


etc.  Next  an  edge  is  taken  from  the  top  of  the  list  and 
it  becomes  the  first  edge  of  the  subgraph.  The  edges  are 
then  sequential ly  taken  from  the  tap  and  added  to  the 
graph  as  long  as  planarity  is  not  violated.  When  3n-6 
edges  have  been  added  the  subgraph  construction  is 
completed.  It  is  noted  that  this  method  requires  an 
explicit  planarity  test  CHopcroft  &  Tarjan,  19741. 

2-2.4  The  N-Boundary  Greedy  Heuristic 

The  N-boundary  Greedy  Heuristic  CGiffin  ft  Foulds, 
19861  is  an  extension  of  the  Greedy  Heuristic  that 
includes  benefits  to  the  final  score  for  not  only 
facilities  that  are  immediately  adjacent  to  one  another 
but  for  facilities  that  are  k  facilities  apart  from  each 
other.  In  addition  to  the  normal  adjacency  matrix 
required,  additional  matrices  that  give  values  for  having 
two  facilities  2,  3,  4,  etc.  facilities  apart  are 
required.  Under  the  assumption  of  approximately  equal 
areas,  normally  a  score  is  higher  if  a  facility  is  fewer 
facilities  distant.  Due  to  this  fact  when  adding  a 
facility  the  shortest  path  to  reach  all  other  facilities 
must  be  calculated  in  order  to  find  an  appropriate 
addition . 
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2-E.5  An  Oriented  Graph  Theoretic  Heuristic 

A  paper  bu  Roth,  Hashimshong,  and  Uachman  C19B23 
suggests  a  method  for  turning  a  graph  into  a  rectangular 
floor  plan,  again  requiring  the  development  of  a  planar 
adjacency  graph.  The  adjacencies  have  na  degree  of 
desirability  in  this  method,  only  a  requirement  for  their 
presence  or  absence.  These  incidence  requirements  are 
converted  into  a  planar  graph  by  the  subtraction  of  edges 
or  the  addition  of  dummy  vertices.  This  planar  graph  is 
then  split  into  tuio  subgraphs  representing  north  south 
and  east  west  orientations  by  a  coloring  technique  and 
dimensions  are  calculated  using  the  PERT  algorithm.  From 
this  technique,  several  alternative  plans  are  generated 
for  further  evaluation.  A  requirement  for  the  dimension 
calculations  is  the  orientation  of  certain  facilities  to 
given  directions.  These  calculations  use  the  PERT 
algorithm  to  find  the  critical  path  from  the  north  side 
of  the  building  to  the  south  as  well  as  a  critical  path 
from  the  west  to  the  east  and  thereby  determine  the 
necessary  dimensions . 

Sr 3  _Qs 1 tahQdcQD -Baaed-Qe tboda 
The  graph  theoretic  heuristics  above  have  a  major 
disadvantage  compared  to  the  Deltahedron  based  heuristics 
that  folloui.  This  disadvantage  is  that  as  yet  there  is 
no  systematic  method  for  finding  the  rectangular 


geometric  dual  to  the  adjacency  graphs  generated.  The 
main  purpose  of  this  thesis  is  to  describe  such  a 
systematic  approach  for  the  deltahedron  based  heuristics. 

A  feature  that  all  of  the  deltahedron  approaches 
have  in  common  is  that  they  begin  with  an  initial 
tetrahedron.  Short  descriptions  of  the  deltahedron 
approaches  follow. 

2-3.1  The  Deltahedron  Heuristic 

The  Deltahedron  Heuristic  CFoulds  and  Robinson, 
197B1  sequentially  adds  a  vertex  into  the  triangle  of  the 
existing  graph  that  will  give  the  greatest  increase  in 
adjacency  scare.  This  increase  in  score  is  calculated  by 
summing  the  weights  of  the  three  edges  used  to  connect 
the  new  vertex  to  the  existing  graph.  The  order  that  the 
vertices  are  added  is  the  continuation  of  the  column  sum 
ordering  used  in  the  initial  tetrahedron  selection  or  the 
selection  at  each  iteration,  of  the  vertex  and  triangle 
that  will  yield  the  greatest  increase  in  score  among  all 
choices  [sometimes  referred  to  as  the  greedy  order!. 

This  method  is  described  in  greater  detail  in  chapter  3 
since  it  is  used  to  generate  the  adjacency  graphs  used  to 
demonstrate  the  development  of  a  block  plan  from  a 
Deltahedron  based  method. 


5-3.2  The  Improved  Deltahedron  Heuristic 

The  Improved  Deltahedron  Heuristic  CFoulds  and 
Robinson,  19793  uses  the  solution  obtained  with  the 
Deltahedron  Heuristic  as  input.  This  graph  is  examined 
to  see  if  any  improvements  can  be  made,  in  the  form  of 
edge  swapping  or  vertex  relocation.  In  most  cases,  if  an 
edge  is  deleted  from  the  graph,  a  quadrilateral  is 
formed.  The  edge  that  was  removed  formed  a  diagonal  in 
this  quadrilateral.  If  the  edge  that  is  identified  with 
the  other  diagonal  is  added  a  new  graph  is  formed  that  is 
also  maximally  planar.  If  the  score  is  increased  by  this 
swap,  it  is  performed,  and  if  not,  it  is  ignored.  All 
possibilities  are  examined  and  when  no  improvements  can 
be  made,  the  process  stops.  In  some  specific  instances 
after  an  edge  is  removed,  the  one  that  would  be  added  is 
already  a  part  of  the  graph.  These  situations  are  either 
ignored,  or  a  well  defined  sequence  of  equivalent  swaps 
made . 

2-3.3  The  N-Boundary  Deltahedron  Heuristic 

As  the  N-Boundary  Greedy  Heuristic  is  an 
extension  of  the  Greedy  Heuristic,  so  too  is  the 
N-Boundary  Deltahedron  Heuristic  CGiffin  ft  Foulds,  19863 
the  same  type  of  extension  to  the  Deltahedron  Heuristic. 
An  increase  to  the  score  of  the  N-Boundary  Deltahedron  is 
determined  by  the  adjacencies  of  facilities  2,  3,  4,  etc. 
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facilities  distant  in  addition  to  the  immediate 
adjacencies  scared  in  the  Deltahedron  Heuristic.  This 
heuristic  begins  with  the  same  initial  tetrahedron 
selection  method  as  the  Deltahedron  method  and  adds  to  it 
by  choosing  the  vertex  that  will  yield  the  highest 
increase  in  score  for  adjacency  or  near  adjacency  to  all 
other  facilities.  As  with  the  N-Boundary  Greedy 
Heuristic,  an  update  version  of  a  shortest  path  algorithm 
must  be  run  at  every  iteration. 

2-3.4  The  Super  Deltahedron  Heuristic 

The  Super  Deltahedron  Heuristic  CGiffin  ft  Foulds, 
19B5D  is  fundamentally  different  from  the  other  graph 
theoretic  methods  in  that  its  objective  function  is  not 
the  maximization  of  total  adjacency  scores;  instead  it 
attempts  to  minimize  transportation  costs  much  like  the 
QAP  formulation  or  the  CRAFT  method.  The  method  again 
starts  with  the  initial  tetrahedron  selection  process 
used  in  the  Deltahedron  method  since  maximizing  the 
proximity  of  four  facilities  with  high  mutual  flows 
should  provide  reasonably  low  transportation  cost.  The 
order  of  insertion  is  either  the  column  sum  or  the  greedy 
approach  used  in  the  Deltahedron  method.  The  triangle 
selected  for  insertion  is  the  one  that  minimizies  the  sum 
of  the  product  of  the  cost  per  unit  distance  traveled, 
the  number  of  trips  per  time  period,  and  the  distance 


between  two  facilities,  over  all  pairs  of  facilities 
contained  in  the  adjacency  graph.  The  shortest  path 
routine  is  also  required  in  this  method  for  the 
computation  of  pairwisB  facility  distances.  ThB  distance 
traveled  between  two  facilities  x  and  y  is  approximated 
by  the  sum  of  half  the  square  root  of  the  area  of  x,  the 
sum  of  the  square  root  of  the  area  of  all  facilities  on 
the  shortest  path  from  x  to  y,  and  half  the  square  root 
of  the  area  of  y.  This  metric  assumes  that  all  facilities 
are  squares  with  side  length  equal  to  the  square  root  of 
the  area,  the  travel  between  two  facilities  is  between 
centroids  of  the  two  facilities,  and  that  all  travel  is 
done  in  a  rectilinear  fashion.  These  assumptions  are  not 
very  likely  in  the  final  block  plan;  however,  they  are 
only  designed  to  give  a  ranking  among  triangles  for  the 
insertion  process. 


CHAPTER  3 


RECTANGULAR  GEOMETRIC  DUAL  AND 
BLOCK  PLAN  CONSTRUCTION 

3=l-lBCBiQQlQgU*-Uo£o£iCHU-QDd_DBfiDifciQD9 
The  Following  terminology  and  notation  is 

def ined : 

Lll-UfiCtBS*  A  point  on  the  adjacency  graph  at 
which  edges  converge  is  known  as  a  vertex. 

ISi.Edgfi •  An  edge  is  a  line  connecting  two 
vertices  on  the  adjacency  graph. 

L32_lnSBCtiQO_QcdBE .  The  insertion  order  is  the 
order  in  which  the  vertices  are  added  to  the  initial 
tetrahedron  to  Form  the  completed  adjacency  graph. 

L22-EBCtQOgulaC-SBQOBfcciC-QuQl-  A  rectangular 
spacial  realization  oF  vertices  and  their  adjacencies 
represented  in  the  adjacency  graph. 

L52_ttQdB»  Each  node  is  a  point  in  the  dual  which 
has  a  one-to-one  correspondence  with  a  triangle  Formed  by 
three  vertices  and  three  edges  oF  the  adjacency  graph. 

L62_Uall>  ^  wall  is  a  line  that  connects  two 
nodes  in  the  dual.  Each  wall  has  a  one-to-one 
correspondence  with  an  edge  in  the  adjacency  graph. 

LZ2_ElQCSd_iD.  When  a  Facility  i  is  added  to  the 

dual,  a  portion  oF  the  dual  is  renamed  to  represent  i. 

SB 
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The  designation  being  replaced  is  called  the  facility 
that  i  uas  placed  in.  If  another  facility  j  was  added  so 
that  a  portion  of  facility  i  is  renamed,  facility  J  is 
placed  in  i,  not  the  original  facility. 

IB2_C0CDfir •  The  right  angle  sometimes  required 
to  connect  two  nodes  of  the  dual  in  a  rectangular  fashion 
is  called  a  corner. 

L92_fidditiQn_Seauence .  ThB  addition  sequence  is 
identical  to  the  insertion  order,  however  it  refers  to 
additions  to  the  dual  not  the  adjacency  graph. 

LlQ2_Uade_ExBansiQD •  Node  expansion  is  the 
redesignation  of  the  structure  surrounding  a  node  in  the 
dual  when  a  facility  is  added  at  that  node. 

Llll_lnhibitac.  An  inhibitor  is  a  dummy  node 
added  to  the  dual  matrix  to  prevent  thB  Idss  of 
adjacencies  when  areas  are  later  added  to  the  dual  to 
form  the  black  plan. 

L1S2_M«  N  is  the  number  of  facilities  or 
vertices . 

L132_ii«.i*.k;>  •  The  combination  of  symbols  <i,J,k> 
represent  the  trianglB  formed  by  vertices  i,  j,  and  k 
with  edges  ij ,  Jk,  and  ki . 

3z£_DeltabedcQD_aetbQd_U3ed 

The  DeltahBdron  Heuristic  seeks  to  find  a 
maximally  weighted  maximally  planar  adjacency  subgraph  of 
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a  complete  adjacency  graph.  The  method  used  here  is  the 
simplest  of  the  variants  of  the  Deltahedrcn  Heuristic. 

The  First  step  is  to  construct  the  NxN  matrix  of  RAj 
values.  The  scores  for  each  Rjj  are  entered  in  the 
matrix  .  The  columns  are  then  summed  and  reordered  in 
nonincreasing  order  by  these  column  sums  with  the 
exception  of  facility  1  which  is  always  the  exterior. 

For  ease  of  discussion,  suppose  that  the  vertices  were 
initially  in  nonincreasing  order  of  column  sums  and 
therefore  their  order  is  1,  2,  ...  ,N.  This  is  now 
refered  to  as  the  Insertion  Order.  ThB  first  four 
verticies  are  combined  to  form  the  complete  graph  on  four 
verticies  K4  which  comprises  the  Initial  Tetrahedron  Csee 
figure  3.1Ca33.  The  vertices  are  then  added  according  to 
the  insertion  order.  Consider  the  insertion  of  vertex  r 
into  triangle  <i,J,k>.  The  benefit  to  the  total  score 
is  the  sum  of  Wir  ♦  Wjr  +  .  Therefore  r  is  chosen  to 

maximize  this  sum  over  all  available  triangles.  After 
adding  vertex  r  into  triangle  <i,J,k>,  this  triangle 
<i,j,k>  is  replaced  by  triangles  <i,j,r>,  <i,k,r>,  and 
<J,k,r>.  The  next  vertex  is  then  selected  and  inserted 
into  the  triangle  that  will  achieve  the  greatest  benefit. 

If  there  is  a  tie  among  several  triangles  for  this 


maximum  benefit,  several  different  strategies  can  be 
incorporated .  One  such  strategy  is  arrange  them 
lexicographically  and  chose  the  first  one.  Another  is  to 
chose  one  of  the  possible  triangles  randomly  and  this 
approach  is  taken  herB  to  avoid  a  large  concentration  of 
insertions  in  one  section  of  the  graph. 

3r3-0esccistiOD-Q£-tbB-BBc£gDaulac-GBOOBtcic.Dual 

Construction 

The  method  used  for  constructing  the  rectangular 
geometric  dual,  hereafter  referred  to  as  the  dual,  is 
limited  to  the  class  of  adjacency  graphs  that  can  be 
constructed  using  any  variant  of  the  deltahedron 
heuristic.  The  only  input  required  is  the  triangle 
insertion  order.  The  process  begins  with  a  rectangular 
representation  of  the  dual  corresponding  to  thB  initial 
tetrahedron.  This  is  shown  in  figure  3.1.  The 
facilities  are  numbered  as  shown  with  facility  1  being 
defined  as  the  exterior.  It  should  be  noted  that  each 
node  of  the  rectangular  geometric  dual  has  three  and  only 
three  edges  incident  with  it.  Each  node  has  a  one  to  one 
correspondence  with  a  triangle  that  exists  in  the 
deltahedron  at  the  current  stage  of  the  adjacency  graph 
construction.  If  a  facility  is  added  to  the  rectangular 
geometric  dual  by  expanding  about  one  of  these  nodes,  its 
adjacencies  will  correspond  exactly  to  those  in  the 
adjacency  graph. 
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< 1 , 2 , 3> 


<  1 , 2^3>Jk<l ,  2 , 4* 
■"  <2,3, 4>^^*>*> 


< 1 . 3, 4> 


<1,3.4> 


Adjacency  Graph 
Cal 


Rectangular  Geometric  Dual 
Cbl 


Figure  3.1.  Initial  Tetrahedron 


There  are  tuo  mays  that  a  facility  may  bB  added 


to  the  dual  with  the  decision  being  made  by  inspection  of 


the  nodes  in  the  dual  that  are  only  one  edge  distant.  If 


there  are  no  corners  that  are  between  thB  node  of 


interest  and  any  of  the  three  adjacent  nodes,  then  the 


facility  is  added  by  a  BOX  operation.  If  there  is  a 


corner  immediately  adjacent  to  the  node  of  interest,  a 


CARUE  operation  is  used.  An  example  of  each  follows. 


3-3.1  Boxing 


From  inspection  of  the  initial  block  plan,  figure 


S.ICbl,  it  can  be  seen  that  the  only  node  that  has  no 


corners  adjacent  to  it  is  <2,3,4>  therefore  consider  the 


insertion  of  facility  5  at  this  node.  From  the  adjacency 
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graph  in  figure  3.2Cal,  it  can  be  seen  that  when  facility 

5  is  inserted  into  <2,3,4>,  the  triangle  <2,3,4>  is 
replaced  by  three  triangles,  namely  <2,3,5>,  <2,4,5>,  and 

<3,4, 5> .  Figure  3.2Cbl  illustrates  this  insertion  and 
the  necessary  relabeling. 
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<2,3,S> 

<2. 4 , 5> 

<1,2,4> 

5 

<3,4,5> 

3 

4 

Cal 

Adjacency  Graph 


<1,3,4> 

CbD 

Rectangular  Geometric  Dual 


Figure  3.2.  Insertion  of  Facility  5  into  triangle 
<2, 3, 4>  CBOXD 


Since  facility  5  replaced  a  portion  of  facility  3,  this 
is  defined  as  placing  facility  5  in  facility  3.  This 
operation  is  called  a  ’’box"  for  obvious  reasons.  The  box 

could  also  be  flipped  to  the  opposite  side  of  the  wall 
separating  facilities  3  and  4.  The  choice  is  arbitrary, 
however  it  does  affect  the  orientation  of  the  block  plan 
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from  the  decision  point  on.  For  ong  given  location  this 
flipping  is  not  aliuags  possible  for  other  reasons  that 
will  be  shown  later.  Four  orientations  of  the  boxing 
operation  arB  possible  and  for  implementation  purposes 
are  defined  as  left-dawn,  left-up,  right-down,  and 
right-up  Csee  figure  3.33. 


Left-Down 


Right-Down 


_n|-  jd 

Left-Up 

jZZLor  JZL 

Right-Up 


Figure  3.3.  Possible  Boxing  Alternatives 


3-3.2  Carving 

Now  consider  instead,  the  insertion  of  facility  5 
into  triangle  <1,2,3>.  This  could  be  done  as  a  boxing 
operation  Cright  and  down]  however  this  would 
unnecessarily  create  an  ”L"  shape  which  is  not  as 
desirable  as  a  rectangle.  This  is  avoided  by  an 
operation  called  a  "carve.”  Figure  3.4  is  an 
illustration  of  this  operation. 


1 


3  4  <  1 , 3, 5>  <  1 , 3, 4> 


Cal  Cbl 

Adjacency  Graph  Rectangular  Geometric  Dual 

Figure  3.4.  Insertion  of  facility  5  into  triangle 
< 1 , 2 , 3>  CCARUE1 


The  same  general  triangle  replacement  is  done  as  above. 
The  eight  orientations  for  the  carve  operation  are  shown 
in  figure  3.5  along  with  their  designations.  These 
designations  indicate  first  the  direction  in  which  the 
corner  is  encountered  followed  by  the  direction  not  cut 
off  by  the  corner.  A  carve  operation  cannot  be  flipped 
to  the  opposite  side  of  a  wall  like  the  box  since  there 
is  no  corner  to  "carve”  towards.  Boxing  might  be  an 
alternative;  however,  as  will  be  shown  later  there  could 


be  a  problem  with  maintaining  the  required  adjacencies  in 
the  dual  when  areas  are  introduced. 


..n.'.'aV.' 


Down-RigKt  Up-Right  Left-Down 


Down-Left 


Up-Left 


Right-Up 


Figure  3.5.  Possible  Carving  Alternatives 


Right-Down 


Left-Up 


Using  these  two  operations,  the  entire  dual  is 
constructed  by  adding  each  facility  to  the  existing  dual 
using  the  same  sequence  used  when  inserting  the  triangles 
in  the  adjacency  graph.  After  the  dual  is  completed,  the 
block  plan  is  made  by  incorporating  the  areas  of  the 
individual  facilities  into  the  orientation  developed 
during  the  dual  construction. 

3z2  _Data  _Stcuc  t  uco  _and  _Camau  tec  _  Lain  lemontatiQD_zz 

DELiaei.au 

The  computer  program  for  this  method  is  called 
DELTAPLAN  and  was  written  in  BA5ICA  on  an  IBM  Personal 


Computer.  Due  to  the  amount  of  memory  available  in 
BA5ICA,  the  problem  size  is  somewhat  limited  however;  11 
facility  problems  can  be  handled  routinely  and  in  some 
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cases  it  will  run  completely  with  as  many  as  S3 
facilities . 


3-4.1  Initialization 

To  facilitate  an  easily  envisioned  and 
manipulated  representation  of  the  dual,  a  matrix  of 
alphanumeric  strings  is  generated  that  contains  the 
elements  common  to  all  initial  block  plans.  As  can  be 
seen  in  figure  3.E,  all  of  thB  initial  triangles  are 
represented  as  six  character  strings.  For  example 
triangle  <1,3,3>  is  represented  by  010E03.  The  walls  are 
represented  by  a  single  dash  ”-”  and  the  interior  of  a 
facility  by  a  two  digit  numerical  string  for  example  ”04” 
for  facility  4  and  ”13”  for  facility  IE.  Since  each 
corner  is  adjacent  to  only  two  facilities  thB  first  two 
elements  of  the  string  are  letters  that  represent  the 
orientation  of  the  corner  Csee  figure  3.7.3  The  two 
corners  in  facility  E  Cupper  left  and  upper  right 
corners]  are  not  used  as  no  facilities  are  added  within 
facility  S  and  therefore  are  represented  by  ’’OOOIOS". 
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Figure  3.6.  Matrix  representation  of  the  rectangular 
geometric  dual 


L_  _i  n  r 

AA  BB  CC  DO 


Figure  3.7.  Corner  Labels 


Ulith  the  exception  of  two,  all  of  the  elements  listed 
above  have  a  direct  counterpart  in  the  dual  shoun  in 
Figure  S.lCbl.  These  exceptions  are  called  "inhibitors” 
and  their  purpose  will  be  defined  later. 


3-4. 8  Addition  of  Facilities  to  the  Rectangular  Geometric 
Dual 

Upon  completion  of  the  deltahedron  heuristic,  for 
simplicity  oil  facilities  are  relabelled  according  to 
their  position  in  the  insertion  order.  Hence,  we  assume 
facility  Ci+43  is  added  to  the  dual  at  the  Cilth  stage 
and  that  facilities  1  through  4  make  up  the  initial 
tetrahedron.  As  can  be  seen  from  figure  3.1Ca3,  only 
four  options  exist  for  the  placement  of  this  first 
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facilitu  and  the  output  of  the  deltahedron  heuristic  used 
to  generate  the  insertion  order  has  chosen  the 
appropriate  one.  A  search  is  then  made  to  match  the 
triangle  in  which  facility  5  is  to  be  inserted,  with  its 
identical  element  in  the  dual  matrix.  A  sort  routine  is 
included  in  the  program  to  insure  consistent  ordering  of 
the  three  two  digit  pairs  within  each  element.  Since  a 
search  of  the  whale  matrix  is  rather  time-consuming,  a 
table  is  constructed  which  contains  each  possible 
insertion  triangle  along  with  its  coordinates  CI,J3  in 
the  matrix. 

£12_SeaccbiDB •  Before  the  search  is  done,  all 
flags  Cdescribed  below!  and  all  direction  indicators  are 
set  to  zero.  Starting  at  the  coordinates  CI,JD,  a  search 
is  performed  to  the  left  to  identify  the  structure  of  the 
dual  to  the  left  of  the  triangle  in  question.  A  variable 
”L”  is  used  to  keep  track  of  the  search  and  is  initially 
equal  to  J .  L  is  decrimented  by  one  and  the  element  with 
coordinates  Cl, LI  is  examined.  If  L  is  less  than  l,  the 
border  of  the  matrix  has  been  reached  and  the  left 
direction  is  "unusable."  An  unusable  direction  means 
that  no  box  or  carve  operation  is  possible  in  this 
direction.  In  the  program  this  is  accomplished  by 
setting  LFLAG0"1 .  If  the  element  is  a  dash,  the 

search  continues  with  the  next  element.  If  a  six  digit 


element  is  found,  the  search  stops.  If  the  first  digit 
of  the  element  is  ”A”  or  ”D”  Cthese  ore  the  only  possible 
corners  when  searching  to  the  left!,  LFLAG1-1  or  4 
respectively.  This  flag  indicates  uihether  a  box  or  a 
carve  operation  is  appropriate  inhere  a  type  A  corner  is 
indicated  with  a  1,  type  B  with  a  2,  type  C  with  a  3,  and 
type  D  with  a  4.  The  presence  of  a  ”000000”  element 
indicates  a  inhibitor  and  the  inhibitor  flag  LFLAG2  is 
set  to  1  C inhibitors  are  described  later  in  this 
chapter.!  If  L-J-l  or  J-2,  the  left  direction  is  again 
unusable  since  there  are  not  enough  elements  between  J 
and  L  to  define  a  new  facility.  After  the  search  to  the 
left,  a  similar  search  is  done  in  the  right,  down,  and  up 
directions. 

LSl-CQCUB^BQX-QBCiSiQD .  The  flags  LFLAG0, 

LFLAG1 ,  RFLAG0,  RFLAG1 ,  DFALG0,  DFLAG1 ,  UFLAG0,  and 
UFLAG1  are  compared  to  the  set  of  values  required  for 
each  carve  operation  to  see  if  it  is  possible  to  carve. 
For  each  carve  operation  three  flags  must  be  set  to 
specific  values.  For  example,  to  carve  left-up  the 
corner  encountered  in  the  left  search  must  be  a  type  A 
CLFLAG1-13,  the  left  direction  must  be  usable  CLFLAG0-0! 
and  the  up  direction  must  be  usable  CUFLAG0-03 .  If  none 
of  the  above  conditions  are  satisfied,  the  flags  required 
for  the  boxing  operations  are  checked.  In  this  case 
there  are  only  two  flags  required  for  a  box  operation. 


For  the  box  left-up  operation  the  flags  needed  are  the 
samB  as  for  the  carve  left-up  CLFLAG0“0  and  UFLAB0-01 
except  there  must  be  no  corner  present  so  the  left  and  up 
corner  indicators  must  be  0  CLFLAG1-0  and  UFLAG1-0D . 

£2J_CaC¥iD0 <  The  left-up  carving  operation  will 
be  used  here  for  description  purposes.  However,  the  same 
general  format  applies  to  all  eight  carving  operations. 
Consider  the  insertion  of  facility  5  into  triangle 
<1,3,4>.  An  inspection  of  figure  3.6  gives  the  structure 
surrounding  010304  and  indicates  that  a  left-up  carve  is 
appropriate.  The  coordinates  CI.JD  of  010304  are 
determined  and  will  become  the  location  of  one  of  the  new 
nodes  of  facility  5.  In  this  case  L  equals  the  j 
coordinate  af  AA0102,  U  equals  the  i  coordinate  of  020304 
and  both  the  right  and  down  directions  are  unusable. 

Next,  the  coordinates  CI1.J1D  of  the  point  diagonally 
across  the  new  facility  from  CI.JD  are  determined.  If 
the  element  which  determines  U  is  not  an  inhibitor,  II  is 
half  way  between  I  and  U.  If  it  is,  Il-U+1,  since  if  an 
inhibitor  is  present,  the  element  above  has  an  unusable 
down  direction.  A  carve  that  goes  only  half  way  up 
wastes  the  entire  portion  above  the  carve  and  is  then 
lost  to  further  insertions.  However,  if  the  carve  goes 
as  close  as  possible  to  the  node  above,  only  a  few 
elements  in  the  matrix  are  lost.  The  J  coordinate  J1  is 


equal  to  L .  In  order  to  determine  the  orientation  of  the 
facilities  which  border  the  new  one,  threB  more  variables 
are  set.  In  this  case  they  are  LS“"01”,  US"”03”,  and 
R$-”04”,  and  they  are  taken  from  the  matrix  by 
determining  which  facilities  are  to  the  left,  right  and 
above  the  new  facility.  These  three  pairs  along  with  the 
number  of  the  new  facility  CFACSD  are  combined  to  form 
the  four  new  nodes  in  the  matrix.  The  upper  right  node 
is  US  +  RS  +  FACS  C03040S]  with  coordinates  CI1.J3,  while 
the  upper  left  nodB  is  LS  ♦  US  +  FACS  C010305D  at 
CI1.J13.  The  lower  left  element  is  ”AA”  ♦  LS  +  FACS 
C AA0105D  at  CI,J1D  and  finally  the  lower  left  node  is  LS 
+  R$  +  FACS  C01040S3  at  C I , J3 .  The  walls  are  then 
inserted  by  renaming  the  elements  between  each  node  on 
the  perimeter  of  the  new  facility  with  The  interior 

of  the  facility  is  then  filled  in  with  FACS  or  in  our 
case  ”05”.  Two  inhibitors  are  then  added  in  place  of  the 
elements  immediately  above  the  upper  left  and  upper  right 
nodes.  The  purpose  of  these  is  described  later.  Figure 
3.B  shows  the  matrix  with  facility  5  added  at  <1,3,4>. 
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Figure  3.8. 
at  < 1 , 3, 4> 
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5 

Tuio  additional  items  are  required  for  the  area 
calculations  that  begin  following  the  completion  of  the 
dual .  The  first  of  these  is  the  operation  with  which  the 
facility  was  added.  In  the  above  example,  the  operation 
is  carve  left  up  therefore  the  variable  OPERSC53 
Coperatian  for  facility  53  is  designated  ”CLU”.  The 
other  requirement  far  the  area  calculations  is  the  number 
of  the  facility  in  which  the  new  facility  was  placed. 

The  variable  for  this  is  PLINS,  and  its  value  in  the 
above  example  is  3  since  the  05  elements  replaced  03 
elements . 

L42_BQSiDg .  The  box  operation  is  accomplished  in 
much  the  same  manner  as  the  carve.  For  this  description, 
the  addition  of  facility  5  at  <2,3,4>  will  be  used.  The 
surrounding  structure  here  indicates  that  a  box  left  down 
operation  is  appropriate.  Notice  that  without  the 
Inhibitor  to  the  right  of  020304  a  box  right  down  would 


We 
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also  be  possibility.  As  noted  earlier,  this  topic  will 
be  discussed  later.  As  in  the  carve  operation,  the 
coordinates  CI.JD  of  020304  are  determined,  as  well  as  L 
and  0,  in  this  case,  L  is  the  J  coordinate  of  010203  and 
D  is  the  i  coordinate  of  010304.  Since  neither  of  these 
is  an  inhibitor,  II  is  half  way  between  I  and  D  and  J1  is 
half  way  between  J  and  L.  If  the  node  to  the  left  had 
been  an  inhibitor,  J1  would  have  been  L+l  and  if  the  node 
below  was  and  inhibitor,  II  would  have  been  D-l .  The 
same  matrix  conservation  reasoning  applies  here  as  in  the 
carve  operation.  The  variables  LS,  US,  and  RS  are  set  as 
described  above  in  order  to  define  the  new  nodes.  Here 
L$-”03”,  US* ”02” ,  RS-”04",  and  FACS  is  again  ”05".  The 
new  nodes  are  020305  for  the  upper  left,  020405  for  the 
upper  right,  030405  for  the  lower  right,  and  AA0305  for 
the  lower  left  element.  As  before,  the  walls  are 
inserted,  interior  of  the  new  facility  is  relabelled, 
0PERSC5D  is  set  to  its  value  of  BLU,  and  PLINSC5D  is  SBt 
to  its  appropriate  value  which  is  03.  A  representation 
of  this  is  given  in  figure  3.S.  It  is  noted  here  that  as 
above  there  are  two  inhibitors,  one  to  the  left  of  the 
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Figure  3.9.  Matrix  representation  with  facility  5  added 
at  <2, 3,4> 

I5J_lDblblt0C8 •  The  purpose  of  inhibitors  is  to 
block  the  Insertion  of  facilities  at  certain  locations 
that  could  possibly  destroy  an  existing  adjacency  once 
areas  are  added.  Consider  the  addition  of  facility  5  to 
<2, 3, 4>  and  the  subsequent  addition  of  facility  6  to 
<2,4, 5> .  If  facility  5  were  added  as  described  above,  it 
is  noted  that  the  coordinates  of  020405  are  the  same  as 
were  the  coordinates  of  020304.  With  the  inhibitor 
present,  as  is  shown  in  figure  3.9,  the  only  possible 
operation  is  a  box  left  down.  However,  if  the  inhibitor 
were  not  present,  a  right  down  box  would  also  be 
possible.  If  the  box  left  down  for  facility  5  were 
followed  by  a  box  right  down  for  facility  6,  the  result 
would  be  as  is  shown  in  figure  3. lOCal.  The  problem 
arises  when  areas  are  introduced.  If  the  area  of 
facility  6  is  larger  than  that  of  facility  5,  the 
adjacency  between  facilities  4  and  5  is  lost  and  an 
adjacency  between  3  and  6  is  gained  as  is  shown  in  figure 
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3.10Cb3.  In  this  case  the  block  plan  would  not  reflect 
the  adjacencies  required  by  the  adjacency  graph.  The 
block  plan  that  does  reflect  the  required  adjacencies 
regardless  of  areas  is  shown  in  figure  3.10Cc3. 

Another  example  of  inhibitors  using  the  carve 
operation  is  illustrated  in  figures  3.10  CdD  and  Cel. 

Here  a  carve  for  facility  5  at  010304  is  followed  by  a 
carve  at  010405  for  facility  6.  With  no  inhibitors,  the 
problem  here  is  the  addition  of  facility  7  at  030405  and 
the  two  options  of  box  left  down  and  box  right  down.  As 
is  seen  in  figure  3.10Cd3  the  box  left  down  destroys  the 
adjacency  between  4  and  5  and  creates  an  adjacency 
between  6  and  7;  however,  at  this  stage  facility  7 
should  only  be  adjacent  to  3,  4,  and  5.  The  box  right 
down  is  appropriate  here  and  figure  3.10Ce3  illustrates 
the  block  plan  which  the  inhibitors  require. 

A  final  example  is  shown  in  figures  3.10  Cf] 
and  CgD.  In  this  case  facility  5  is  added  at  030304 
followed  by  a  carve  for  facility  6  at  030405.  When 
facility  7  is  added  at  030305,  the  same  problem  presented 
in  figure  3.11  arises  again.  With  no  inhibitors  the 
block  plan  could  end  up  as  in  figure  S.lOCfl,  whereas 
inhibitors  require  the  black  plan  in  figure  3.10Cg3. 

The  initial  choice  of  location  for  the 
inhibitors  to  the  right  of  030304  and  010304  is 
arbitrary.  Placement  of  both  on  the  left  would  perform 
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just  as  well  but  it  should  be  noted  that  they  must  both 
be  on  the  same  side  or  thou  would  create  the  very 
problems  they  ore  designed  to  eliminate. 

The  results  below  follow  from  the  operations  as 

def ined . 

L62_IbBQCam_l.  No  more  than  one  carve  can  be 
done  within  any  Facility.  PROOF  —  In  order  to  carve 
there  must  be  a  corner  towards  which  one  carves.  After 
one  carve  is  done,  there  is  no  corner  left  in  the 
original  facility  therefore  the  condition  required  to 
carve  does  not  exist  and  no  further  carving  can  be  done. 

IZ2_IbBDCBiD_£ .  No  more  than  three  facilities  may 
be  placed  within  any  given  facility  i.  PROOF  —  All 
facilities,  with  the  exception  of  8,  begin  as  boxes. 

Even  if  a  facility  is  added  by  a  carve  operation  it 
contains  one  corner  and  therefore  has  the  same  structure 


Figure  3.11.  Location  of  inhibitors  when  no  Facilities 
may  be  added 


IB2_CQEQliQru_2.il  •  IF  three  facilities  are  added 
within  facility  i,  two  must  be  boxes  and  one  a  carve. 
PROOF  —  For  a  given  node,  if  there  is  an  opportunity  to 
carve  it  will  be  done  first.  From  theorem  1,  one  cannot 
carve  again  therefore  the  other  facilities  must  be  added 
by  a  box  operation. 

From  Corollary  2.1,  the  worst  shape  a  facility 
may  have  is  a  ”T” . 

3-4:. 3  Creating  The  Block  Plan 

The  block  plan  is  nothing  more  than  addition  of 
areas  to  the  dual.  To  accomplish  this  it  is  easiest  to 
start  with  a  ’’clean  slate”  rather  than  trying  to  adjust 
the  existing  dual.  The  inputs  required  for  each  facility 
i  in  this  phase  are  the  operation  COPERSCill,  the 
facility  that  it  was  placed  in  CPLINSCill,  and  the  area 
CAREAC ill.  Each  facility  in  the  block  plan  is  given  by 
its  coordinates  within  a  square  with  sides  of  length  one 


and  where  the  coordinates  represent  percentages  of  the 
actual  wall  lengths.  For  example,  consider  two  buildings 
each  containing  10,000  square  feet,  with  dimensions 
100x100  for  the  first  and  125x80  for  the  second  Csee 
figure  3.12.3  ft  facility  with  dimensions  CO, 03, 

CO, 0. S3 , CO. 5,03 ,  and  CO. 5, 0.53  would  have  dimensions  of 
50x50  in  the  first  case  and  62.5x40  in  the  second  however 
as  one  can  see  the  areas  are  both  equal  to  2,500  sq.  ft. 
This  adds  mare  flexibility  to  the  actual  site  block  plan 
since  no  restriction  is  made  that  the  building  be  square. 
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•  125- 
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40 
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Figure  3.12.  Coordinate/ArBa  Relationship 

IlJ_CoiD0utlD0-tbB.lDltlol-6cBO-BBayiced.foi:-Bocb 
Facility-  The  area  required  for  a  facility  i  when  it  is 
initially  added  into  the  block  plan  is  not  the  area  of 
facility  i  alone  since  subsequent  facilities  are  added 
within  the  initial  boundaries  of  facility  i.  The  initial 
facility  should  contain  the  area  required  for  all  of  the 


facilities  added  uithin  its  initial  boundaries  at  later 
stages.  Using  thB  PLIN  vector,  a  cumulative  area  vector 
called  AREAIN  is  calculated  so  that  each  value  of 
AREAINCi3  is  equal  to  the  area  of  facility  i  plus  the 
cumulative  areas  of  all  facilities  subsequently  added 
uithin  the  initial  boundaries  of  facility  i. 

I£J_CDryiDg_ODd_IbB_0dditiDD_Df_tbo_Eirst_IwD 
Eacilities_£Q_£be_BlDCls_ElBD  •  The  entire  square  is 
defined  as  the  initial  boundary  of  facility  5,  therefore 
its  cumulative  area  CAREAINCE33  is  equal  to  the  total 
area  or  AREATOT .  Facility  3  is  then  placed  uithin  the 
initial  boundary  of  facility  S.  Since  the  initial 
facility  3  contains  all  facilities  except  2  it  can  be 
vieued  as  a  carve  up  from  belou.  It  is  noted  that  both 
the  carve  left  up  and  the  carve  right  up  look  the  same 
uith  the  only  difference  being  the  node  from  uhich  the 
carve  took  place.  In  the  initial  dual  section  this  uas 
an  important  distinction,  houever  for  the  block  plan  it 
doesn’t  really  matter  since  the  shape  for  the  block  plan 
is  all  ue  are  concerned  uith  here  Csee  figure  3.5.3 
Therefore  in  the  block  plan  section  only  four  carve 
routines  are  required  since  the  left-up  and  right-up,  the 
left-doun  and  right-doun,  the  doun-right  and  up-right, 
and  the  doun-left  and  up-left  are  equivalent.  The  carve 
operation  at  this  stage  involves  bosically  cutting  the 
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initial  area  of  2  into  tuo  parts  that  have  the  proper 
ratio  of  areas.  Since  the  coordinates  are  in  percentages 
of  distance,  the  carve  operation  mag  be  accomplished  bg 
simply  relabeling  the  louer  coordinates  of  facilitg  2  as 
the  louer  coordinates  of  the  initial  area  of  facilitg  3, 
redefining  the  louer  tuo  coordinates  of  facilitg  2 
according  to  the  ratio  of  cumulative  areas,  and  also 
assigning  these  coordinates  as  the  upper  coordinates  of 
the  initial  area  of  facilitg  3.  The  cumulative  area  of 
facilitg  3  CAREAINC333  is  then  subtracted  from  the 
cumulative  area  of  2  CAREAINC233  to  get  the  neu 
cumulative  area  of  facilitg  2.  The  same  tgpe  of 
operation  is  done  for  adding  the  initial  area  of  facilitg 
4  uithin  facilitg  3  but  a  carve  to  the  to  left  is  used. 

Up  to  this  paint  there  have  been  no  problem 
specific  facilities  placed  as  facilities  2  through  4 
aluags  have  the  same  initial  location.  From  here  on,  the 
facilities  are  not  necessarilg  added  in  the  same 
sequence  as  theg  uere  in  the  insertion  order;  instead 
theg  are  added  according  to  the  facilitg  that  theg  are 
placed  in.  For  example,  all  facilities  uhose  PLIN  value 
is  3  are  added  to  facilitg  3,  then  those  uith  PLIN  values 
of  4,  etc.  From  Theorem  2  and  its  Corollarg,  at  most 
three  facilities  mag  be  placed  in  facilitg  i  and  theg 
must  be  a  subset  of  tuo  boxes  and  a  carve.  The  PLIN 
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exist,  that  are  placed  in  facility  i.  If  a  carve 
operation  is  present,  it  is  done  first.  The  carve  method 
described  above  for  the  initialization  of  facilities  2 
through  4  is  used  for  subsequent  carve  additions. 

I3i_BQS_fiddi£iQDS_£Q_£bB_BlQCk_El0D .  When  there 
are  two  boxes  to  be  added  to  the  block  plan,  the  one  with 
the  largest  cumulative  area  is  chosen  to  be  inserted 
first.  Consider  the  addition  of  facility  5  at  <E,3,4> 
within  facility  3  as  described  above  Csee  figure  3.2.3 
The  upper  right  coordinates  of  facility  3  are  relabeled 
as  the  upper  right  coordinates  of  facility  5.  The  lower 
right  and  upper  left  coordinates  are  calculated  according 
to  the  square  root  of  the  ratio  of  cumulative  areas.  The 
lower  left  coordinate  is  the  i  coordinate  of  the  lower 
right  and  the  j  coordinate  of  the  upper  left.  The  only 
change  to  the  existing  facility  C33  is  relabeling  of  the 
upper  right  coordinate  which  is  the  same  as  the  lower 
left  of  the  new  facility. 

The  addition  of  a  second  box  is  done  in  the  same 
manner  as  the  first  so  long  as  there  is  sufficient  space. 
If  there  is  not  a  "correction”  routine  is  entered.  The 
definition  of  "sufficient  space”  is  as  follows.  After 
one  box  has  been  added,  an  L  shape  exists.  The 
coordinates  for  the  rectangular  portion  of  this  existing 
L  shape  where  the  new  box  is  to  be  added  are  used  to 


determine  the  "effective”  area  of  the  existing  facility. 
If  the  area  of  the  box  to  be  added  is  more  than  95%  of 
this  effective  area,  there  is  not  sufficient  space.  If 
this  is  the  case,  wall  length  of  the  first  box  in  the 
offending  direction  is  reduced  with  the  adjacent  wall 
being  increased  to  maintain  the  specified  area.  When 
sufficient  space  is  achieved,  the  second  box  is  added 
along  with  the  corrected  first  box.  As  with  the  dual 
construction,  these  operations  are  used  repeatedly  until 
the  block  plan  is  completed. 


CHAPTER  4 

EXAMPLE  PROBLEMS 

In  this  chapter  DELTAPLAN  solutions  to  three 
different  problems  are  presented.  The  first  example  is  a 
problem  from  Francis  and  White  C1974D  and  comparisons 
uiith  ALDEP  and  CORELAP  solutions  are  given.  The  second 
example  is  also  from  Francis  and  White,  and  it  includes 
the  illustration  of  a  possible  extension  to  include 
changes  to  the  adjacency  graph  made  bg  the  Improved 
Deltahedron  Heuristic.  The  final  example  is  a  problem 
that  is  too  large  to  be  solved  bg  the  current  version  of 
DELTAPLAN,  houiever  a  brief  description  of  the  variable 
reassignment  required  to  construct  the  complete  block 
plan  is  included. 

4=l-£xQ(DDle-l 

The  first  example  is  a  ten  facility  problem 
houiever,  since  the  Deltahedron  method  requires  the 
exterior  to  be  included  as  facility  1,  the  problem  shown 
has  11  facilities.  The  REL  chart  required  as  input  by 
the  Deltahedron  method  is  given  in  figure  4.1. 
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_ Eccility^ezes _ 

1  Exterior  - V. 


E  Display  Area 

- ISO  -  — < 

3  Parcel  post 

- 340  -  -- < 


4  Parts  shipment 


5  Foreman 


6  Repair  and  service  parts 


7  Service  area 


B  Receiving 


9  Testing 


10  General  storage 


11  General  offices 


Figure  4.1.  Example  I  REL  Chart 


The  insertion  order  calculated  using  column  sums  is: 


1  10B724956  11  3 


From  the  insertion  order  it  can  be  seen  that  the  initial 
tetrahedron  is  1-10-0-7  and  table  4.1  gives  the  remaining 
vertices  and  the  triangles  into  which  they  were  inserted. 


Table  4.1  Example  I  Vertices  and  Insertion  Triangles 

_ Vectex _ Iciangle _ 

2  <  1  B  7> 

4  <  1  10  7> 

9  <10  B  7> 

5  <10  7  9> 

6  <  5  7  9> 

11  <2  B  7> 

3  <10  7  5> 


Figure  4.3.  Example  I  Block  Plan 

The  complete  actual  output  from  this  example  is 
given  in  the  appendix.  In  addition  to  the  output  given 
here,  the  appendix  includes  the  incidence  matrix,  a 
condensed  version  of  the  AS  matrix,  the  insertion  order 
information,  and  the  coordinates  of  the  block  plan.  The 
incidence  matrix  is  a  duplicate  of  the  original  REL  chart 
with  the  adjacencies  not  present  in  the  adjacency  graph 
replaced  by  dashes.  The  condensed  AS  matrix  uses  numbers 


59 


to  represent  the  interior  of  facilities,  dashes  to 
represent  the  walls  C including  intersections],  and  0’s  to 
represent  the  inhibitors.  The  first  line  of  the 
insertion  ordBr  information  gives  thB  second,  third,  and 
fourth  facilities  inserted,  and  their  areas.  Each 
additional  line  gives  the  facility  number,  the  area,  the 
operation  used  to  insert  the  facility  in  the  dual,  the 
triangle  it  was  placed  in  Crelabeled  to  correspond  to  the 
order  of  insertion],  and  the  facility  that  the  new 
facility  was  placed  in  Calso  relabeled].  The  coordinates 
listed  are  in  the  same  relative  position  on  the  page  as 
in  the  block  plan  i.e.  the  upper  left  coordinate  of  each 
group  of  four  is  the  coordinate  of  the  upper  left  corner 
of  the  facility.  In  the  case  where  a  box  has  been  placed 
in  a  facility  and  there  are  now  six  corners  in  the 
facility,  the  coordinate  of  the  corner  where  the  box  was 
placed  is  the  coordinate  of  the  box  that  protrudes  into 
the  old  facility  Csee  figure  4.4]. 


Figure  4.4.  Coordinate  location  when  a  Box  is  placed 
within  the  facility 

Figures  4.5  and  4.5  show  the  output  From  ALDEP 
and  CORELAP  for  the  same  problem.  For  comparison,  the 
scares  for  each  are  calculated  using  the  scoring  rules  of 
the  Deltahedran  method.  This  is  Justified  since  the 
scoring  for  the  ALDEP  method  is  identical  Cthis  is  true 
in  this  case  since  there  are  no  facilities  adjacent 
diagonal lyl  and  CORELAP  includes  maximization  of 
adjacencies  in  its  objective  function.  5cores  for 
adjacencies  to  the  exterior  are  not  included  since  the 
input  for  ALDEP  and  CORELAP  solutions  did  not  include 
these  adjacencies  in  their  REL  charts,  therefore  the 
scores  for  adjacency  with  the  exterior  are  subtracted  off 
the  Oeltahedron  score. 
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Figure  4.5.  ALDEP  Layout  for  Example  I 
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Comparison  shows  that  the  Deltahedron  method 
achieved  the  highest  score  with  217  followed  by  ALDEP 
with  211  and  Finally  CORELAP  with  210.  It  is  noted  that 
there  are  several  narrow  L  shaped  facilities  in  the 
DELTAPLAN  block  plan  however  some  modifications  described 
in  the  next  example  and  in  chapter  5  might  help  to  create 
more  rectangular  or  regular  spaces. 

2=S_EsaarlB_il 

The  second  example  is  also  from  Francis  and  White 
and  as  with  example  I,  an  exterior  facility  has  been 
added  resulting  in  a  12  facility  problem.  The  PEL  chart 
for  this  problem  has  been  rearranged  so  that  the 
insertion  order  is  simply  increasing  integers  from  1  to 
12.  Figure  4.7  gives  the  REL  chart  used  for  input  and 
table  4.2  gives  the  insertion  vertices  and  triangles. 
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Since  the  REL  chart  has  been  rearranged,  the  initial 
tetrahdron  is  1-2-3-4. 


_ EOCllitu/eCBQ _ SstlDB 


Figure  4.7.  Example  II  REL  Chart 


Table  4.2  Example  II  Vertices  and  Insertion  Triangles 
_ UeCtBJS _ ICiQDBlO _ 

5  <1  2  4> 

6  <2  3  4> 

7  <5  2  4> 

6  <2  3  6> 

9  <2  3  B> 

10  <2  3  9> 

11  <1  2  3> 

12  <1  4  S> 


Figure  4.9.  Example  II  Block  Plan 

Consider  the  general  manager’s  office  Cfacilitg 
51  and  note  that  the  L  shape  is  not  desirable. 

Consulting  the  REL  chart  it  is  also  noted  that  the 
adjacency  rating  between  facilities  7  and  IE  is  an  E  and 
the  adjacency  rating  between  7  and  4  is  only  a  U.  With  a 
series  of  edge  swaps  of  the  type  described  in  the 
improved  deltahedron,  an  increase  in  score  can  be 
achieved  while  also  making  the  general  manager’s  office  a 


iginal  insertion  order 
CaD 


e 


u 


Adjacency  graph  and  dual  after  three  edge  swaps 

Cd3 

Figure  4.10 — ContiDUBtl. 

It  should  be  noted  that  since  there  is  no  vertex  in  the 
subgraph  illustrated  with  degree  three,  there  is  no 
possible  way  to  generate  this  graph  using  only  the 
Deltahedron  method  as  the  last  vertex  inserted  must  have 
degree  three.  Additionally,  the  current  Improved 
Deltahedron  would  not  consider  this  sequence  of  changes 
since  the  first  swap  results  in  a  lower  score;  therefore 
a  look  ahead  procedure  would  be  required.  It  is 
therefore  proposed  that  every  permissible  edge  swap  can 
be  characterized  in  the  dual  Cand  the  block  plan!  as 
transforming  a  carve  into  a  box  or  a  box  into  a  carve. 

It  is  further  proposed  that  since  every  maximally  planar 
graph  can  be  constructed  from  an  initial  tetrahedron  by  a 
series  of  vertex  insertions  and  edge  swaps,  CGiffin, 

19043  if  the  sequence  is  known  it  is  possible  to 
construct  the  dual  of  all  maximally  planar  adjacency 
graphs.  The  computer  implementation  of  this  procedure 
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has  not  yet  been  done  nor  has  the  multiple  step  look 
ahead  implementation  of  the  Improved  Deltahedron  method. 
This  is  left  For  future  research. 

S-3_ExamDle_III 

The  final  example  is  a  real  world  problem  and 
illustrates  the  degeneracy  that  often  occurs  in  some 
larger  actual  problems.  It  also  illustrates  the  outcome 
of  a  problem  that  is  too  large  to  be  solved  by  the 
current  program.  Consider  the  REL  chart  illustrated  in 
figure  4.11 . 
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It  is  clear  from  an  inspection  of  the  REL  chart  that 
facilities  2,  5,  11,  and  IB  have  no  rating  other  than  □ 
with  respect  to  all  other  facilities.  Therefore  the 
problem  is  degenerate  because  when  choosing  among 
triangles  for  facilities  2,  5,  11,  or  IB,  any  triangle  is 
as  good  as  every  other.  Additionally,  facilities  3  and  4 
have  an  E  only  between  themselves  and  an  0  with  all 
others  as  do  B  and  7.  It  follows  that  as  long  as  3  and  4 
are  adjacent  and  E  and  7  are  adjacent,  a  block  containing 
facilities  2,  3,  4,  5,  6,  7,  11,  and  16  could  be  placed 
anywhere  in  the  graph  and  result  in  the  same  score  as 
placing  it  anywhere  else.  Because  Df  this  property, 
there  are  literally  thousands  of  combinations  that  would 
result  in  the  same  score  but  have  different  adjacency 
graphs.  One  approach  to  this  dilemma  might  be  to  group 
the  B  facilities  into  one  large  facility  and  thus  reduce 
the  size  of  the  problem  by  more  than  a  third.  For  the 
sake  of  demonstration  however,  the  entire  problem  is  run 
as  given.  This  illustrates  the  problem  encountered  by 
the  current  program  when  the  AS  matrix  becomes  too  small 
to  add  all  of  the  required  facilities  within  it.  The 
Deltahedron  method  runs  without  incident  with  the  initial 
tetrahedron  being  facilities  1-19-22-21  and  the  insertion 
vertices  and  triangles  are  given  in  table  4.3.  As  with 
examples  I  and  II,  the  complete  output  is  contained  in 
the  appendix. 
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Table  4.3  Example  III  Uertices  and  Insertion  Triangles 
- ysrtsjs _ IriODole _ 


B 

<19 

EE 

Sl> 

IE 

<  1 

19 

31> 

10 

<  B 

EE 

El> 

9 

<19 

El 

B> 

13 

<  1 

EE 

El> 

19 

<19 

EE 

B> 

EO 

<  1 

19 

1£> 

3 

<  8 

El 

10> 

4 

<  B 

El 

3> 

5 

<  B 

10 

3> 

7 

<  5 

10 

3> 

14 

<E0 

19 

1E> 

15 

<1B 

EE 

B> 

S 

<  B 

10 

B> 

5 

<  B 

5 

E> 

11 

<  1 

El 

13> 

16 

<  B 

3 

4> 

17 

<15 

EE 

B> 

The  resultant  dual  is  given  in  Figure  4. IE  and  it  should 
be  noted  that  there  are  only  19  Facilities  shown. 
Facilities  E,  5,  and  15  were  not  able  to  be  inserted 
since  there  was  no  room  at  the  new  location  For  an 
additional  Facility.  The  program  can  be  continued 
normally  From  this  point  and  output  obtained,  however  the 
block  plan  will  not  contain  thB  Facilities  leFt  out  oF 
the  dual  Csee  Figure  4.133. 


Figure  4.13.  Example  III  Block  Plan  with  3  facilities 
not  included 

To  provide  a  complete  block  plan,  the  BREAK  feature  of 
BASIGA  is  used.  Before  continuing  on  to  the  construction 
of  the  block  plan  from  the  dual,  the  program  execution  is 
stopped  with  the  BREAK  keg .  When  the  program  is  halted 
in  this  manner,  the  variables  defined  up  to  this  point 
remain  in  memory.  The  values  not  present  for  the 
complete  construction  of  the  block  plan  are  the  variables 
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□PERS  and  PLIN  for  facilities  2,  5,  and  16.  An 
inspection  of  the  condensed  AS  matrix  along  with  the 
insertion  values  displayed  on  the  screen  yield  the 
necessary  information  to  determine  mhat  the  values  mould 
have  been  had  the  program  had  the  necessary  room.  In 
this  case  the  folloming  variables  mere  set  to  the  values 
indicated  belarn. 


□PERSC1B1 
□PERSE  191 
□PERS C 21 D 


BRD” 

PLINC1B3-14 

CDL” 

PLINC191-1B 

CDR” 

PLINC211-13 

After  these  values  are  set,  execution  is  resumed  and  the 
result  is  given  in  figure  4.14.  It  should  be  noted  that 
for  different  random  seed  values,  DELTAPLAN  mill  complete 
this  problem  mith  no  variable  redefinition  required. 


2  5  6  4 

16 


Figure  4.14.  Example  III  Block  Plan  complete 


A  possible  change  to  allow  somewhat  larger 
problems  would  be  to  redefine  the  data  structure  and  to 
have  an  AS  matrix  that  starts  out  very  small  and  expands 
out  from  the  point  where  an  additional  facility  is 
placed.  This  is  opposite  to  the  present  method  which 
starts  with  a  given  size  and  facilities  are  placed  within 
its  boundaries.  Many  times  there  is  quite  a  lot  of  space 
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CHAPTER  5 


CONCLUSIONS 

The  method  presented  here  has  been  shouin  to 
construct  □  rectangular  geometric  dual  for  the  class  of 
adjacency  graphs  developed  using  Deltahedron  based 
heuristics.  It  has  further  demonstrated  hou  areas  can  be 
incorporated  to  form  a  block  plan.  It  should  be 
remembered  houiever,  that  all  of  the  methods  described  in 
this  thesis  are  analytical  in  nature  and  as  stated  in 
Francis  and  White  C19743,  "It  should  be  realized  that  the 
analytical  approach  yields  a  solution  to  the  model,  but 
not  necessarily  the  pcgblBID.”  For  this  reason,  one 
should  be  cautious  uhen  selecting  a  block  plan  produced 
by  any  of  the  heuristics  mentioned.  Just  because  a 
particular  plan  has  a  higher  adjacency  score  does  not 
mean  that  it  is  a  better  plan.  The  maximally  planar 
plans  developed  by  DELTAPLAN  sometimes  have  long  narroui  L 
or  T  shaped  facilities  which  are  most  likely  not  very 
useful  if  included  as  shown  in  the  block  plan.  The 
output  of  this  as  well  as  other  methods  is  meant  to  be  a 
starting  place  and  guide  for  further  planning. 
Alternatives  that  may  not  previously  have  been  considered 
might  surface  with  a  computer  method  such  as  DELTAPLAN. 


1 


|'V,V,V 


As  a  starting  point,  the  plan  and  REL  chart  mag  be 


consulted  to  see  if  perhaps  one  of  the  adjacencies  in  a 


long  narrow  L  or  T  shape  is  even  worthwhile  to  have  and 


as  such  it  may  become  a  candidate  for  deletion.  If  a 


graph  is  not  very  dense  in  highly  weighted  edges,  perhaps 


a  maximally  planar  block  plan  contains  more  adjacencies 


than  are  really  necessary.  In  this  case  some 


adjacencies  may  bB  deleted  to  form  a  more  regular  plan 


and  the  adjacency  score  may  not  even  be  affected. 


A  very  important  fact  is  that  the  dual  is  not 


unique.  There  are  many  ways  of  arranging  facilities  with 


very  slight  changes  to  the  rules  of  DELTAPLAN,  that 


preserve  all  of  the  adjacencies  required.  One  change 


might  involve  moving  the  initial  inhibitors  from  the 


right  side  of  the  wall  between  facilities  3  and  4  to  the 


left  side.  Another  possibility  would  be  to  change  the 


placement  of  facilities  2,  3,  and  4  within  the  dual 


representation  of  thB  initial  tetrahedron  which  would 


lead  to  six  different  orientations  of  the  initial  four 


facilities.  These  are  either  2,  3,  or  4  on  the  top  and 


the  remaining  two  facilities  placed  either  on  the  left 


and  right  or  the  right  and  left.  The  point  here  is  that 


the  rules  developed  for  DELTAPLAN  continue  to  work  for 


all  of  these  orientations.  If  areas  are  not  a  factor  or 


if  it  can  be  determined  that  no  new  facility  might  affect 


certain  adjacencies  after  areas  are  added,  changes  to  the 


wivivu 
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inhibitors  at  a  later  stage  can  be  invoked  an  result  in 
still  further  alternatives. 

As  mentioned  earlier,  extensions  can  be  made  to 
include  edge  swaps  of  the  type  used  in  the  Improved 
Oeltahedron.  Further  extensions  include  the  ability  to 
develop  the  dual  from  any  maximally  planar  adjacency 
graph  once  the  series  of  vertex  insertions  and  edge  swaps 
required  to  form  the  adjacency  graph  from  an  initial 
tetrahedron  are  known .  As  yet  it  is  not  known  how  the 
process  of  enforcing  a  deltahedron  like  insertion  and 
swapping  procedure  on  an  arbitrary  adjacency  graph  should 
be  performed  efficiently.  Another  extension  might  be  to 
develop  the  block  plan  in  parallel  with  the  Super 
Oeltahedron  method  in  order  to  have  more  accurate 
estimates  of  the  distance  between  centroids  for 
transportation  cost  estimates. 

The  implementation  provided  in  this  thesis  should 
form  an  important  subroutine  to  the  realization  of  all  of 
these  extensions. 


E 


APPENDIX  A 


THE  DELTAHEDRDN  HEURISTIC 
PROGRAM  LISTING 


10 

ao 

30 

40 

50 

60 

70 

80 

90 

100 

110 

iao 

130 

140 

150 

160 

170 


180 

190 

aoo 

eio 

aao 

if  not 

330 

340 

350 

360 

370 

380 

390 

300 

310 

330 

330 

340 

350 

360 


The  DELTAHEDRON  HEURISTIC 
using  column  sums  insertion  order 

by  J.  W.  Biffin 

with  modifications  by  D.  Ul.  Keenan 
flarch  1,  1986 


Din  BEN  C30. 303. ORDER  C303 , BENSUM  C303 , TRIANG  C64 , 33 , SOLUTION  C30.303 
Din  OTHERS  C303.DEG  C303.A  C303 , DEBCON  C303.R0DTA  C303 
Din  SPATH  C 30, 303 , HP  C33.H  C303 , UALIDC303 , TRI ANGSC303 , 

AREAC303 , PSC30, 303 
CLS 

RANDOMIZE 
DEFINT  I-N 

INPUT  "You  uiill  need  to  input  the  Filename  for  the  data  you  want  to 
use.  Would  you  like  a  list  of  files  on  the  disk 

CY/N3 ” ; ANSS 

IF  ANSS- "Y"  OR  ANSS-"y"  THEN  FILES 

INPUT“Enter  any  filename  with  .DAT  For  an  extension” ; FILENAMES 
IF  RIGHTSCFILENAMES.430~.DAT”  THEN  F I LENAnES-F I LENAMES* ” . DAT " 

PRINT 

INPUT” If  you  need  an  X  value  other  than  -1034  enter  it  at  the  prompt, 
press  return . ”  ; XUALS 

IF  XUALS-”"  THEN  XUAL— 1034  ELSE  XUAL-UALCXUALS3 

I  _ _ _ _  _  _  _  _ 

'Read  data  from  data  file  and  initialize  NxN  score  matrix 


OPEN  ” I " , #1 , FILENAMES 
INPUT  *1,  N 

PRINT  "NUMBER  Or  FACILITIES: ”N 

FLAB  -0 

K-7 

FOR  1-1  TO  N 

PRINT  USING  "»#” ; I ; : PRINT  ”  :  " 
PRINT  TABCK3 

PRINT  USING  ”**";!  ;:PRINT  ”  "j 
FOR  J-I+l  TO  N 


B1 


rn>. 

VS* 

lT« 


370 

300 


450 

460 

470 

400 

490 

500 

510 

520 

530 

540 

550 

560 

570 

500 

590 

600 

610 

620 

630 

640 

650 

660 

670 

600 

690 

700 

710 

720 

730 

740 

750 

760 

770 

700 

790 

000 

810 


INPUT  *1.  PSCI.JD 
PRINT  PSC  I  , J  1  " 


390 

IF  PSCI, J 3-”U”  THEN 
450 

BEN 

CI.J3-0  : 

BEN  CJ.ID-0  : 

GOTD 

400 

IF  PSC I , J3-"0”  THEN 
450 

BEN 

CI.J3-1  : 

BEN  CJ.I3-1  : 

GOTO 

410 

IF  PSCI,  J3-’'I”  THEN 
450 

BEN 

CI.J3-4  : 

BEN  C J, 13-4  : 

GOTO 

420 

IF  PSCI . J3--E”  THEN 
450 

BEN 

CI.J3-16  : 

BEN  CJ, 13-16 

:GOTO 

430 

IF  PSC  I  ,J3-”A’‘  THEN 
450 

BEN 

C I , J  3-64  : 

BEN  C J , I 3-64 

:GOTO 

440 

IF  PSCI, J3-”X"  THEN 

BEN 

C I , J3-XUAL 

.  :  BEN  CJ.I3- 

XUAL 

NEXT  J 
K-K+2 
PRINT 
NEXT  I 
FOR  1-2 

NEXT  I 
CLOSE 


:  FLAG 


TO  N 

INPUT  #1,  AREACI3 


'If  an  X  is  present,  add  a  constant  to  all  scores  so  thou 
'are  all  non-negative 

—  ______  —  ______ ,  _  _  __ _ _ _  _  _ _ 

IF  FLAG  -0  GOTO  630 
FOR  1-1  TO  N 

FOR  J-l  TO  N 

BEN  C  I ,  J  3-BEN  CI.J3-XUAL 

NEXT  J 

NEXT  I 

FOR  1-1  TO  N 

BEN  Cl, 13-0 

NEXT  I 
GOSUB  1530 

'Initialize  total  score  and  add  the  score 
'For  the  initial  tetrahedron 

TOTBEN  -0 
FOR  1-1  TO  4 

FOR  J-I+1  TO  4 

TOTBEN  -TOTBEN  *BEN  C ORDER  C I  3 , ORDER  CJ33 

NEXT  J 

NEXT  I 

’Determine  best  triangle  for  vertex  insertion 

FOR  1-5  TO  N 

MAX— 1 


020 

830 

040 

050 

060 

070 

080 

090 

900 

910 

920 

930 

940 

950 

960 

970 

980 

990 

1000 

1010 

1020 

1030 

1040 

1050 

1060 

1070 

1080 

1090 

1100 

1110 

1120 

1130 

1140 

1150 

1160 

1170 

1100 

1190 

1200 

1210 

1220 

1230 

1240 

1250 

1260 


X-OROER  cn 

CX-1 ♦ I NTC  RND-TR 1 NO  1 

FDR  X-CX  TO  TR1N0 

sun  -0 

FOR  J-l  TO  3 

sun  -sun  +ben  cx  .triang  cx,J33 

NEXT  J 

if  sun  >  n«x  then  nAX-sun  :  iiaxtri  -k 

NEXT  X 

FOR  X— 1  TO  CK-1 
sun  -o 

FOR  J-l  TO  3 

sun  -sun  +ben  cx  .triang  cx.jn 

NEXT  J 

if  sun  >  nftx  then  nAX-sun  :  nAXTRi  -x 

NEXT  X 


‘Print  vertex  and  triangle  it  is  inserted  into 


PRINT  “INSERTING  VERTEX  ”;X  IN  TRIANGLE 
FOR  X-  1  TO  I 

IF  TRIANG  C ttAXTRI , 1 3  -  ORDER  CX3  THEN  ELttNTl-X:  GOTO 
1050 

NEXT  X 

PRINT  TRIANG  C ttAXTRI  ,13; 

FOR  X-  1  TO  I 

IF  TRIANG  C ttAXTRI , 23  -  ORDER  CX3  THEN  ELNNT2-X :  GOTO 
1090 

NEXT  X 

PRINT  TRIANG  C nAXTRi  ,23; 

FOR  X-  1  TO  I 

IF  TRIANG  C ttAXTRI , 33  -  DRDER  CX3  THEN  ELCINT3-X :  GOTO 
1130 

NEXT  X 

PRINT  TRIANG  C ttAXTRI  ,33; 


'Create  character  sting  elements  used  as  input  for  DELTAPLAN 


IF  ELttNTl< 10  THEN  ELttNTlS-"0”+RIGHTSCSTRSCELnNTl 3,13  ELSE 
ELttNT 1 S-R I GHTS  C  STRS  C  ELttNT 13,23 

IF  ELHNT2< 10  THEN  ELttNT2S-"0"*RIGHTSCSTRSCELttNT23 , 1 3  ELSE 
ELttNT2S-RIGHTSC5TRSCELttNT23 ,23 

IF  ELHNT3<10  THEN  ELttNT3S-’‘0”*RIGHTSCSTRSCELttNT33 , 1 3  ELSE 
ELttNT 3S-RIGHTSCSTRSCELMNT33, 23 
TRIANGSC I 3-ELttNTlS*ELttNT2S+ELttNT3$ 

PRINT  TRIANGSC 1 3 
PRINT 

GOSUB  2130 
TOTBEN  -TOTBEN  -ttAX 

NEXT  I 

IF  FLAG-1  THEN  TOTBEN-TOTBEN  «■  XVAL»C3»N-63 
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1270 

1280 

1290 

1300 

1310 

1320 

1330 

1340 

1350 

1360 

1370 

1380 

1390 

1400 

1410 

1420 

1430 

1440 

1450 

1460 

1470 

1480 

1490 

1500 

1510 

1520 

1530 

1540 

1550 

1560 

1570 

1580 

1590 

1600 

1610 

1620 

1630 

1640 

1650 

1660 

1670 

1680 

1690 

1700 

1710 

1720 

1730 

1740 

1750 


PRINT 

PRINT  "TOTAL  DELTAHEDRON  ADJACENCY  SCORE  IS”  TOTBEN 
PRINT 

GOSUB  2280 

'Writs  output  to  data  File 

• _ .  _  _  _  - _ .  -  ..  . _ --  _ 

OPEN  ”0" , #1 , "DATA1 ” 

WRITE  #1 , N 
FOR  1-1  TO  N 

WRITE  ttl , ORDERC I } 

NEXT  I 

FOR  1-5  TO  N 

WRITE  ttl.TRIANGSIIl 

NEXT  I 

FOR  1-2  TO  N 

WRITE  #1, AREA C ORDERC II 3 

NEXT  I 
CLOSE 

INPUT  "WOULD  YOU  LIKE  A  LAYOUT  DONE  FROn  THIS  DATA  CY/N3”;ANSS 
IF  ANSS- "N"  OR  ANSS-"n”  GOTO  1490 
CHAIN  "DELTAPLN" 

END 


Print  NxN  scorn  matrix 


FOR  1-1  TO  N 

FOR  J-l  TO  N 

PRINT  BEN  CI.J3; 

NEXT  J 
PRINT 

NEXT  I 
PRINT 


Calculate  column  sums 


FOR  J-l  TO  N 

SUM  -0 
FOR  1-1  TO  N 

IF  IOJ  THEN  Sun  -SUn  +BEN  CI.J3 

NEXT  I 

BENSun  CJi-sun 

NEXT  J 

FOR  1-1  TO  N 

UAL ID  Cl  3-1 

NEXT  I 

FOR  I-  1  TO  N 


ORDER  Cl  3-1 


NEXT  I 


B4 


1760 

1770  '  Sort  vertices  according  to  column  sums 

17B0  '  Bubblesort  array  order  according  to  BENSUtl 

1790 

1B00  FLIPS  -1 

1810  WHILE  FLIPS  -1 

1880  FLIPS  -0 

1830  FOR  I ”2  TO  N-l 

1840  IF  BENSUn  CORDER  C 1 3  3  <  BENSUfl  CORDER  C 1+111  THEN  SWOP 

ORDER  C I  3 , ORDER  CI  +  13  :FLIPS  -1 

1850  NEXT  I 

1860  WEND 

1870 

1880  'Print  deltahedran  insertion  order 

1890 

1900  PRINT 

1910  PRINT  "DELTAHEDRON  INSERTION  ORDER” 

1980  PRINT 

1930  FOR  1-1  TO  N 

1940  PRINT  ORDER  CI3; 

1950  NEXT  I 

1960  PR I NT : PR I NT 

1970 

1980  'Initialize  triangles  and  incidence  values  for  the 

1990  'initial  tetrahedron 

8000 

8010  FOR  1-1  TO  4 

8080  X  -ORDER  CI3 

8030  FOR  J-l  TO  4 

8040  Y  -ORDER  CJ3 

8050  IF  J<I  THEN  TRIANG  CI.J3-Y  ELSE  IF  J> I  THEN  TR1ANG 

CI.J-13-Y  :  SOLUTION  CX , Y 3-1 : SOLUTION  CY.X3-1 

8060  NEXT  J 

8070  NEXT  I 

8080  TRINO  -4 

8090  RETURN 

8100 

8110  '<<<  Relabel  deleted  triangle  and  odd  tuo  more  >>> 

8180 

8130  FOR  J-l  TO  3 

8140  SOLUTION  CX  .TRIANG  CflAXTRI ,  J3 3-1 

8150  SOLUTION  CTRIANG  CnAXTRI,J3,X  3-1 

8160  NEXT  J 

8170  TRINO  -TRINO  ♦! 

8180  TRIANG  CTRINO  ,13-TRIANG  CflAXTRI, 13 

8190  TRIANG  CTRINO  ,83-TRIANG  CflAXTRI, E3 

8800  TRIANG  CTRINO  ,33-X 

8810  TRINO  -TRINO  +1 

8880  TRIANG  CTRINO  ,13-TRIANG  CflAXTRI  ,  1 3 

8830  TRIANG  CTRINO  ,83-TRIANG  CflAXTRI, 33 

8840  TRIANG  CTRINO  ,33-X 


1 1  > 


m 


B5 


5550 

5560 

5570 

5500 

5590 

5300 

5310 

5350 

5330 

5340 

5350 

5360 

5370 


TRIANG  CHAXTRI.ll-X 
RETURN 


•<<<  Print  matrix  of  adjacencies  present  >>> 


PRINT  "INCIDENCE  flATR I X :  “ 

PRINT 

K-7 

FOR  1-1  TO  N 

PRINT  I; 

PRINT  TABCK3 
FOR  J-I+l  TO  N 

IF  SOLUTION  CI.Jl-l  THEN  PRINT  PSCI.JJ; 


ELSE  PRINT 


5300 

5390 

5400 

5410 

5450 

5430 


NEXT  J 

K-K*S 

PRINT 


NEXT  I 

PRINT 

RETURN 
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10 

so 

30 

40 

50 

60 

70 

BO 

90 

100 

110 

ISO 

130 

140 

150 

1E0 

170 

180 

190 

SOO 

S10 

SSO 

S30 

S40 

SSO 

seo 

S70 

SBO 

S90 

300 

310 

3S0 

330 

340 

350 

360 

370 

3B0 


DELTAPLAN 

A  Procedure  to  Construct  a  Rectangular  Geometric  Dual  and 
a  Block  Plan  from  a  Deltahedron  Based  Adjacency  Graph 

by  David  W.  Keenan 
(larch  1,  1986 


KEY  OFF 
SCREEN  1 
WIDTH  BO 
CLS 


— tl.JJ  ARE  THE  INSERTION  TRIANGLE  COORDINATES 

--111.J1J  ARE  THE  OPPOSITE  CORNER  COORDINATES  OF  THE  NEW  FACILITY 
— L  IS  THE  HORIZONTAL  COORDINATE  OF  THE  NEXT  ADJACENCY  TRIANGLE  TO 
THE  LEFT 

— R  IS  THE  HORIZONTAL  COORDINATE  DF  THE  NEXT  ADJACENCY  TRIANGLE  TO 
THE  RIGHT 

— D  15  THE  UERTICAL  COORDINATE  OF  THE  NEXT  ADJACENCY  TRIANGLE  BELOW 
— U  IS  THE  UERTICAL  COORDINATE  OF  THE  NEXT  ADJACENCY  TRIANGLE  ABOUE 
— _FLAG0-0  INDICATES  THE  DIRECTION  IS  USABLE 
— _FLAG0-1  INDICATES  THE  DIRECTION  IS  NOT  USABLE 

— _FLAG1  INDICATES  A  CORNER  OR  THE  ABSENCE  THEREOF  —  O-NO  CORNER, 
1-LOWER  LEFT  CORNER,  S-LOWER  RIGHT  CORNER,  3-UPPER  RIGHT  CORNER, 
4-UPPER  LEFT  CORNER 


DIH  AS C 35, 66} , RISC SOO 1 , RSCS003 , R3C5001 , PLINC503 , PLINSI503 , 0PERSC50J , 
TR I ANGLES  C  50  3 , AREA  C  50  3 , AREA I NC  50  3 , ORDER  C  503 
PLACE  -  3 
CT-4 

BLS-STR I NGSC 79,333 
’>>>  INITIALIZE  hATRIX 
GOSUB  1040 

•>>>  GET  INPUT  AND  DETERHINE  ITS  COORDINATES 
FOR  FAC-5  TO  NUttFAC 
L-0 
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400  D-0 

410  U-0 

420  1-0 

430  J-0 

440  IF  FAC  <  10  THEN  FACS-“0"+RIGHTSCSTRSCFAC3 . 1 3  ELSE 

FACS-R I GHTS  C  SIRS  C FAC 1,53 

450  GOSUB  5430  ’>>>  CLEAR  FIRST  LINE 

460  PRINT  "Inserting  Facility” ;0RDERCFAC1 ; 

470  LOCATE  2,1 

480  PRINT  BLS 

490  AREA INC FAC 3-AREA C FAC 3 

SOO  ’>>>  SORT  THE  INPUT  TO  INSURE  PROPER  CHARACTER  SEQUENCE 

S10  SORTAS-TR I ANGLES  C  FAC  3  :  GOSUB  5070  :  TR I ANGLES  C  F AC  3 -SORT AS 

550  FOR  X-l  TO  CT 

530  IF  TR I ANGLES! FAC 3 -RISC X 3  THEN  I-RSCX3 : J-R3CX3  :R2CXl-0 

:R3CX3-0  : GOTO  550 

540  NEXT  X 

550  IF  1-0  THEN  LOCATE  5,1  :  PRINT  “THIS  TRIANGLE  CAN’T  BE  FOUND 

AS  LISTED  —  TRY  AGAIN  ”:GOTO  BBO 
560  ’>>>  BEGIN  SEARCH 

570  •>>>  SEARCH  LEFT 

5B0  GOSUB  1870 

590  •>>>  SEARCH  RIGHT 

600  GOSUB  5010 

610  ’>>>  SEARCH  DOWN 

620  GOSUB  5150 

630  ’>>>  SEARCH  UP 

640  GOSUB  5290 

650  ’>>>  CHECK  FOR  CORNERS  AND  CARUE  IF  PQ5SIBLE 


660 

IF  LFLAG1-1 
820 

AND 

LFLAGO-O 

AND 

UFLAGO-O 

THEN 

GOSUB 

5430 

:G0T0 

670 

IF  LFLAG1-4 
820 

AND 

LFLAGO-O 

AND 

DFLAGO-O 

THEN 

GOSUB 

5640 

:  GOTO 

680 

IF  RFLAG1-B 
BBO 

AND 

RFLABO-O 

AND 

UFLAGO-O 

THEN 

GOSUB 

2B50 

:  GOTO 

690 

IF  RFLAG1-3 
820 

AND 

RFLAGO-O 

AND 

DFLAGO-O 

THEN 

GOSUB 

3060 

:GOTO 

700 

IF  DFLAG1-1 
850 

AND 

DFLAGO-O 

AND 

RFLAGO-O 

THEN 

GOSUB 

3570 

:  GOTO 

710 

IF  DFLAG1-B 
850 

AND 

a 

r 

2> 

G) 

O 

1 

o 

AND 

LFLAGO-O 

THEN 

GOSUB 

34B0 

:  GOTO 

750 

IF  UFLAG1-3 
BBO 

AND 

UFLAGO-O 

AND 

LFLAGO-O 

THEN 

GOSUB 

3690 

:  GOTO 

730 

IF  UFLAG1-4 

AND 

UFLAGO-O 

AND 

RFLAGO-O 

THEN 

GOSUB 

3900 

:GOTO 

B50 

740  ’>>>  CHECK  FOR  INTERSECTIONS  AND  BOX  IF  POSSIBLE 

750  IF  LFLAGO-O  AND  DFLAGO-O  THEN  GOSUB  4110  :GOTO  BBO 

760  IF  LFLAGO-O  AND  UFLAGO-O  THEN  GOSUB  4350  :GOTO  BBO 

770  IF  RFLABO-O  AND  DFLAGO-O  THEN  GOSUB  4590  :GOTD  BBO 

7B0  IF  RFLAGO-O  AND  UFLAGO-O  THEN  GOSUB  4B30  :GOTO  BBO 

790  LOCATE  5,1 


800 


PRINT  “This  triangle  cannot  be  inserted  here,  try  onother 
location" 

810  GOTO  880 

820  ’CONTINUE 

830  LOCATE  PLACE, SO 

840  PLINCFAC3-UALCPLINSCFACD3 

850  PRINT  0RDERCFAC3 ; AREACFAC3 ; "  ” ; 0PERSCFAC3 ; "  " ; TRIANGLESCFAC1 ; 

” ; PLINCFAC3 

BBO  IF  PLACE  >-23  THEN  PLACE  -  4  ELSE  PLACE  -  PLACE  «-l 

B70  LINE  C J*4*2B , I*4*503-C J1*4+2B , I 1 •4*503 , , B 

BBO  NEXT  FAC 

890  ’FOR  U-l  TO  CT 

900  1  PRINT  R1SCU3;R2CU3;R3CU3 

910  ’NEXT  U 

920  ’PRINT  CT 

930  GOSUB  5430  ’>>>  CLEAR  FIRST  LINE 

940  INPUT  ’’Would  you  like  a  layout  copy  printed  CY/N3”;AN5S 

950  IF  ANSS  -  "Y”  OR  ANSS  -  “y”  THEN  GOSUB  1770 

9B0  GOSUB  5430  ’>>>  CLEAR  FIRST  LINE 

970  INPUT  "Would  you  like  an  insertion  order  copy  printed  IY/N3’’;ANSS 
980  IF  AN5S  -  "Y"  OR  ANSS  -  “y"  THEN  GOSUB  5370 

990  GOSUB  5430  ’>>>  CLEAR  FIRST  LINE 

1000  INPUT  "Would  you  like  to  see  the  layout  uith  areas  CY/N1”;AN5S 

1010  IF  ANSS--Y”  OR  ANSS-"y"  THEN  GOSUB  5480 

1020  CHAIN  "DELTASUn” 

1030  END  ’ - 

1040  ’>>>  INITIALIZE  OATRIX 

1050  OPEN  ” I " . *1 , "DATA1 " 

10B0  INPUT  #1 ,  NUflFAC 

1070  FOR  1-1  TO  NUMF AC 

1080  INPUT  ei.OROERCU 

1090  NEXT  I 

1100  FOR  1-5  TO  NUMF AC 

1110  INPUT  #1, TRIANGLESC 13 

1120  NEXT  I 

1130  FOR  1-2  TO  NUMF AC 

1140  INPUT  W1.AREACI3 

1150  NEXT  I 

1160  CLOSE 

1170  LOCATE  10,5 

1180  PRINT  "Please  wait  a  Few  moments  while  things  are  being 

initialized . . . . " 

1190  FOR  J-0  TO  BB 

1200  ASCO, J3-"01 “ 

1210  ASC1.J3-”-" 

1220  ASC3.J3-*-" 

1230  ASC34.J3-"-" 

1240  ASC35, JD-"01" 

1250  NEXT  J 

1260  FOR  1-1  TO  34 
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FOR  Jl-2  TO  42 

ASCI ,J13-"03” 

NEXT  J1 

FOR  J2-44  TO  64 

ASC I , J23-”04” 

NEXT  J2 

NEXT  I 

ASC1,13-"000102" 

ASC 1 , 653-”000102" 

ASC3, 1 3“”010203" 

ASC 3 ,43 3 -”020304” 

AS  C  3 , 65  3 - ”010204 " 

ASC34, 13-"AA0103" 

ASC34 , 433-”010304” 

ASC34 , 653-"BB0104” 

ASC 3. 44 3 -”000000” 

AS  C  34 , 44  3- ”000000 “ 

CLS 

LOCATE  PLACE, 50 

PRINT  OROERC23 ; " : " ; AREAC23 ; ORDERC33 : ” : ” ; AREAC33 ;0RDERC4 3 ;  ' 

PLACE-PLACE* 1 

LINE  C32, 543-C288 , 1B63 , , B 

LINE  C32, B23-C200 , 1B63 , , B 

LINE  C200.623-C2BB, 1863 , , B 

RISC  1 3-”010203" 

R2C13-3 
R3C 13-1 

R1SC2 3- "020304" 

R2C23-3 

R3C23-43 

R1SC33-"0X0204” 

R2C33-3 

R3C33-65 

R1SC4 3- "010304" 

R2C43-34 

R3C43-43 

AREA INC  2  3-AREAC  2  3 
AREAINC33-AREAC33 
AREA I NC4 3-AREAC 4 3 
RETURN 

•>>>  <<<  PRINT  AS  MATRIX 
FOR  J-66  TO  0  STEP  -1 


* ; AREAC4  3 


A'l 


VS 


90 


1790 

FOR  1-0 

TO  35 

1800 

PRINT  ASCI.JD  "  "i 

1B10 

IF  LENCASCI , J33-6  AND  ASC I , J )<> "000000”  THEN  AAS 
ELSE  AAS-RIGHTSCASC I , J] , 21 

1820 

IF  AAS-”-”  THEN  LPRINT  ”  ELSE  LPRINT  USING 

“##" ; ORDERCUALC AAS1 3 ; 

1830 

NEXT  I 

1840 

LPR1NT 

1850  NEXT  J 

1050  RETURN 

1870  *>>>  <<<  SEARCH  LEFT 

1880  LFLAGO-O 

1890  LFLAG1-0 

1900  LFLAG2-0 

1910  L-J 

1920  L-L-l 

1930  IF  L<1  THEN  LFLAG0-1:  GOTO  2000 
1940  AUALS-LEFTSCASCI.Ll.ll 

1950  IF  AUALS-”-"  THEN  GOTO  1920 

1960  IF  AUALS-”A”  THEN  LFLAG1-1 

1970  IF  AUALS-”D"  THEN  LFLAG1-4 

1980  IF  AS C I, LJ- "000000”  THEN  LFLAG2-1 
1990  IF  L-J-l  OR  L-J-2  THEN  LFLAGO-1 

2000  RETURN 

2010  ’»>  <<<  SEARCH  RIGHT 

2020  RFLAGO-O 

2030  RFLAG1-0 

2040  RFLAG2-0 

2050  R-J 

2060  R-R+l 

2070  IF  R>65  THEN  RFLAGO-1:  GOTO  2140 

2080  AUALS-LEFTSC ASC I , R3 , 1 3 

2090  IF  AUALS-"-”  THEN  GOTO  2060 

2100  IF  AUALS- “B”  THEN  RFLAG1-2 

2110  IF  AUALS- "C"  THEN  RFLAG1-3 

2120  IF  AS  C I , R ] - “000000 ”  THEN  RFLAG2-1 
2130  IF  R-J+l  OR  R-J *2  THEN  RFLAGO-1 
2140  RETURN 

2150  '>>>  <<<  SEARCH  DOWN 

2160  DFLAGO-O 

2170  DFLAG1-0 

2180  DFLAG2-0 

2190  D-I 

2200  0-0+1 

2210  IF  D>34  THEN  DFLAGO-1 :  BOTO  2280 

2220  AUALS-LEFTSCASCD, J3 , 13 

2230  IF  AUALS-”-"  THEN  GOTO  2200 

2240  IF  AUALS- "A”  THEN  DFLAG1-1 

2250  IF  AUALS- “B”  THEN  DFLAG1-2 

2260  IF  AS C 0 , J 1 - ”000000 ”  THEN  DFLAG2-1 
2270  IF  O-I+l  OR  D-I +2  THEN  OFLAGO-1 
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8280  RETURN 

8230  •>>>  <<<  SEARCH  UP 

8300  UFLAGO-O 

8310  UFLAG1-0 

8380  UFLAG8-0 

8330  U-I 

8340  U-U-l 

2350  IF  U<1  THEN  UFLAG0-1:  GOTO  8480 
83G0  AUALS-LEFTSCASCU, J3 , 13 

8370  IF  AUALS-”-“  THEN  GOTO  8340 

8300  IF  AUALS-'C"  THEN  UFLAG1-3 

8330  IF  AUALS-”0”  THEN  UFLAG1-4 

8400  IF  AS C  U ,  J  3  - ”000000 ”  THEN  UFLAG8-1 
8410  IF  U-I-l  OR  U-I-8  THEN  UFLAGO-1 
8480  RETURN 

8430  ’>>>  <<<  CARUE  LEFT-UP 

8440  IF  UFLAG8-1  THEN  Il-U+1  ELSE  I 1-I-C INTCABSC  J3/833 

8450  Jl-L 

8460  LS-ASCI1+1, Jl-13 

8470  US-ASCU-1,  Jl+13 

8480  RS-ASCI1+1, J+13 

8430  SORTAS-LS+US+FACS  :  GOSUB  5070  : 

AS  C 1 1 , J 1 3 -SORTAS : CT— CT+1 : R1SCCT3 -SORTAS : R8  C  CT  3 - 1 1 : R3  C  CT  3 - J 1 
8500  SORTAS-RS+US+FACS  :  GOSUB  5070  : 

AS t  U , J 3-SORT AS : CT-CT  *1 : R1SICT  3 -SORTAS : R8 1 CT  3 - 1 1 : R3 1 CT  3 - J 
8510  SORTAS-LS+RS+FACS  :  GOSUB  5070  : 

AS  C I , J  3 -SORT AS : CT-CT+ 1 : R 1 S  C  CT  3 -SORT AS : R8  C  CT  3 - 1 : R3  C  CT  3 - J 
8580  SORTAS- ”AA ”-»LS+FACS  :  GOSUB  5070  :  ASC I , J1 3-SORTAS 
8530  ASC 1 1— 1 , J  3  — "000000” 

8540  ASC 1 1-1, J1 3— '000000” 

8550  0PERSCFAC3--CLU" 

8560  PL I NSC FAC 3 -US 

8570  FOR  J8-J1+1  TO  J-l 

8580  ASCII, J83-"-" 

8530  FOR  18-11+1  TO  1-1 

8600  ASC 18, J83-FACS 

8610  NEXT  18 

8680  NEXT  J8 

8630  RETURN 

8640  •>>>  <<<  CARUE  LEFT-DOUN 

8650  IF  DFLAG8-1  THEN  Il-D-1  ELSE  I 1-I+C INTC ABSC I-D3/83 3 

8660  Jl-L 

8670  LS-ASCI1-1, Jl-13 

8680  OS-ASC I 1+1 , Jl+1 3 

8630  RS-ASC I 1-1 , J+l 3 

2700  SORTAS-LS+DS+FACS  :  GOSUB  5070  : 

ASCII , J1 3-SORTAS : CT-CT+1 : R1SCCT 3-SORTAS: R8CCT3-I 1 : R3CCT3-J1 
8710  SORTAS-RS+DS+FACS  :  GOSUB  5070  : 

ASC 1 1 , J  3 -SORTAS : CT-CT+ 1 : R1 S  C  CT  3 -SORTAS : R8  C  CT  3 - 1 1 : R3  C  CT  3 - J 
8780  SORTAS-LS+RS+FACS  :  G05UB  5070  : 

AS  C I , J  3 -SORTAS : CT-CT+ 1 : R1 S  C  CT  3 -SORTAS : R8  C  CT  3 - 1 : R3  C  CT  3 - J 
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I 

1 


m 


8730 

8740 

8750 

8760 

8770 

8780 

8790 

8800 

8810 

8880 

8830 

8840 

8850 

8860 

8870 

8880 

8890 

8900 

8910 

8980 

8930 

8940 

8950 

8960 

8970 

8980 

8930 

3000 

3010 

3080 

3030 

3040 

3050 

3060 

3070 

3080 

3090 

3100 

3110 

3180 

3130 

3140 

3150 

3160 

3170 


SORTAS-"DD"*LS*FACS  :  GOSUB  5070  :  ASCI , J1 3-SORT AS 
ASCI1+1, J3-”000000” 

ASC I 1+1 , J13--000000* 

OPERSC  FAC3““CLD” 

PL I NSC FAC 1 -OS 
FOR  J8-J1+1  TO  J-l 

ASCII, J83-”-" 

FOR  18-1*1  TO  I 1-1 

ASC 18, J83-FACS 

NEXT  18 

NEXT  J8 
RETURN 

•>>>  <<<  CARUE  RIGHT-UP 

IF  UFLA68-1  THEN  Il-U+1  ELSE  I 1-I-C INTCABSC I-Ul/83 3 
Jl-R 

LS-ASC 11*1, J-l 3 
US— ASC 11-1 , Jl-13 
RS-ASC 11*1 , Jl*l 3 

SORT AS— RS+US* FACS  :  GOSUB  5070  : 

AS  C 1 1 . J 1 3 -SORT AS : CT-CT+1 : R1SCCT3 -SORT AS :R8CCT3-I1:R3CCT3-J1 
50RTAS-LS*US*FACS  :  GOSUB  5070  : 

ASCII, J 3-SORT AS : CT-CT* 1 : RISC  CT  3-SORTAS : R8  C  CT  3 - 1 1 : R3  C  CT  3 - J 
SORT AS-LS+RS+F ACS  :  GOSUB  5070  : 

AS  C  I ,  J  3  -SORT  AS :  CT-CT*  1 :  R 1 S  C  CT  3  -SORT  AS :  R8  C  CT  3  - 1 :  R3  C  CT  3  -  J 
SORTAS-”BB”*RS*FACS  :  GOSUB  5070  :  AS  C I , J 1 3 -SORTAS 
ASC 11-1 , J3-"000000“ 

ASC I 1-1, J13--000000" 

0PERSCFAC3--CRU" 

PL I NSC  FAC 3 -US 
FOR  J8-J+1  TO  Jl-1 

ASCII, J83-"-” 

FOR  18-11*1  TO  1-1 

ASC 18, J83-FACS 

NEXT  18 

NEXT  J8 
RETURN 

■>>>  <<<  CARUE  RIGHT-DOWN 

IF  0FLAG8-1  THEN  U-D-l  ELSE  Il-I+C INTCABSC I-D3/83 3 
Jl-R 

LS-ASC I 1-1, J-l 3 

DS-ASC 1 1  +  1 , Jl-1 3 

RS-ASC I 1-1, J1+1D 

SORT AS-RS+DS+F ACS  :  GOSUB  5070  : 

ASCII, J1 3-SORTAS :CT-CT*1: R1SCCT3-SDRT AS: R8CCT3-I 1 : R3CCT3-J1 
SORT AS-LS*DS*F ACS  :  GOSUB  5070  : 

ASCII , J 3-SORTAS : CT-CT+1 : R1SCCT3-S0RTAS: R8CCT3-I 1 : R3CCT3-J 
SORTAS-LS+RS+FACS  :  GOSUB  5070  : 

ASC I , J 3-SORTAS : CT-CT+1 : R1SCCT3-S0RTAS : RSCCT3-I : R3CCT3-J 
SORTAS-”CC”*RS*FACS  :  GOSUB  5070  :  ASC I . J1 3-SORTAS 
ASC I 1*1 , J3»”000000" 

ASC I 1*1 , J1 3- "000000" 
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3100  OPERS  C  F  AC  3 - ”CRD" 

3190  PL  I NSC  FAC 3 “OS 

3500  FOR  J2-J+1  TO  Jl-1 

3210  ASCI1.J2D  —  -" 

3220  FOR  I2-I+1  TO  11-1 

3230  ASCI2, J23-FACS 

3240  NEXT  12 

3250  NEXT  J2 

3250  RETURN 

3270  ’>>>  <<<  CARUE  DOWN-RIGHT 

3280  1 1-0 

3290  IF  RFLAG2-1  THEN  Jl-R-1  ELSE  Jl-J+C INTCABSC J-R3/23 3 

3300  US-ASC 1-1 , Jl-13 

3310  RS-ASCI1-1, Jl+13 

3320  DS-ASCI1+1, Jl-13 

3330  SORTAS-RS+DS+FACS  :  GOSUB  5070  : 

AS  C  1 1 .  J1 3  -SORTAS :  CT-CT+1 :  R1S  t  CT  3  -SORTAS :  R2  C  CT  3  - 1 1 :  R3  C  CT  3  -  J 1 
3340  SORTAS- MAA"+DS+FACS  :  GOSUB  5070  :  ASC II , J3-S0RTAS 

3350  SORTAS-US+DS+FACS  :  GOSUB  5070  : 

ASCI , J 3-SORT AS: CT-CT+1 : R1SCCT3 -SORTAS :R2CCT3-I:R3CCT3-J 
3360  50RTAS-RS+US+FACS  :  GOSUB  5070  : 

ASC I , J1 3-SORTAS: CT-CT+1 : R1SCCT3-S0RTAS: R2CCT3-I : R3CCT3-J1 
3370  ASCI , Jl+1 3-“000000" 

3300  ASCII, Jl+13— 000000” 

3390  0PERSCFAC3— CDR” 

3400  PL I NSC FAC 3 -RS 

3410  FOR  12- I +1  TO  I 1-1 

3420  AS  CI2.J13— 

3430  FOR  J2-J+1  TO  Jl-1 

3440  ASC 12, J23-FACS 

3450  NEXT  J2 

3460  NEXT  12 

3470  RETURN 

3480  ’>>>  <<<  CARUE  DOWN-LEFT 

3490  1 1-0 

3500  IF  LFLAG2-1  THEN  Jl-L+1  ELSE  Jl-J-C INTCABSC J-L3 /23 3 

3510  US-ASC 1-1, Jl+13 

3520  LS-ASC I 1-1, Jl-13 

3530  DS-AS  CI1+1, Jl+13 

3540  SORTAS-LS+DS+FACS  :  GOSUB  S070  : 

AS  C  1 1 ,  J 1 3  -SORT  AS :  CT-CT*  1 R1 S  C  CT  3  -SORTAS :  R2  C  CT  3  - 1 1 1 R3  C  CT  3  -  J 1 
3550  SORTAS- ”BB"+DS+FACS  :  GOSUB  5070  :  ASC I 1 , J3-S0RTAS 
3560  SORTAS-US+DS+FACS  :  GOSUB  5070 

AS  C I , J  3 -SORTAS : CT-CT+ 1 : R 1 S  C  CT  3 -SORTAS : R2  C  CT  3 - 1 : R3C  CT  3 - J 
3570  SORTAS-LS+US+FACS  :  GOSUB  5070  : 

AS  C I , J 1 3 -SORTAS :CT-CT+1:R1SCCT3-S0RTAS:R2CCT3-I:R3CCT3-J1 
3500  ASCI , Jl-1 3—000000" 

3590  ASCII, Jl-13— 000000” 

3600  OPERSC  FAC3-"CDL” 

3610  PLINSCFAC3-LS 

3620  FOR  I 2-1+1  TO  I 1-1 
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3630  ASCI2.J13-”-" 

3640  FOR  J2-J1*1  TO  J-l 

3650  ASC12, J23-FACS 

3660  NEXT  J2 

3670  NEXT  12 

3660  RETURN 

3690  ’>>>  <<<  CARUE  UP-LEFT 

3700  I 1-U 

3710  IF  LFLAG2-1  THEN  Jl-L+1  ELSE  Jl-J-C INTCABSC J-L3/23 3 

3720  US-ASCI1-1, Jl*13 

3730  LS-ASCI1+1, Jl-13 

3740  OS” AS  CI*1,J1*13 

3750  SORTAS-LS+US+FACS  :  GOSUB  5070  : 

AS  C 1 1 , J 1 3 -SORT AS : CT -CT + 1 : R1 S  C  CT  3 - SORT AS : R2  C  CT  3 - 1 1 : R3  C  CT  3 - J 1 
3760  SORTAS- ”CC”*US*FACS  :  GOSUB  5070  :  AS  C 1 1 , J  3 -SORTAS 
3770  SORTAS-US*OS*FACS  :  GOSUB  5070  : 

AS  C I , J  3 -SORTAS : CT-CT* 1 : R1 S  C  CT  3 -SORTAS :  R2  C  CT  3  - 1 :  R3  C  CT  3  -  J 
3780  SORTAS-LS*DS+FACS  :  GOSUB  5070  : 

AS  C I , J 1 3 -SORTAS : CT-CT ♦  1 : R 1 S  C  CT  3 -SORTAS : R2 1  CT  3 - 1 : R3C  CT  3 - J 1 
3790  ASCI, Jl-1 3— 000000“ 

3800  ASCII, Jl-1 3-”000000” 

3810  OPERS  C  F  AC  3 - “CUL “ 

3820  PL I NSC FAC 3 -LS 

3830  FOR  12-11*1  TO  1-1 

3B40  ASCI2, J13-"-" 

3850  FOR  J2-J1+1  TO  J-l 

3860  ASCI2, J23-FACS 

3870  NEXT  J2 

3880  NEXT  12 

3890  RETURN 

3900  ’>>>  <<<  CARUE  UP-RIGHT 

3910  I 1-U 

3920  IF  RFLAG2-1  THEN  Jl-R-1  ELSE  J1-J*C INTC ABSC J-R3/23 3 
3930  US-ASCI1-1, JX-13 

3940  RS-ASC 11*1 , Jl*l 3 

3950  DS-ASC I*1,J1-13 

3960  SORTAS-RS*US*FACS  :  GOSUB  5070  : 

AS  C 1 1 .  J 1 3  -SORTAS :  CT-CT*  1 :  R1 S  C  CT  3  -SORTAS :  R2  C  CT  3  - 1 1 :  R3  C  CT  3  -  J 1 
3970  SORTAS- "DD”*US*FACS  :  GOSUB  5070  :  ASC I 1 , J3-S0RTAS 

3980  SORTAS-US*DS*FACS  :  GOSUB  5070  : 

AS  C I , J  3 -SORTAS : CT-CT* 1 : R 1 S  C CT  3-SORTAS : R2  C  CT  3 - 1 : R3  C  CT  3 - J 
3990  SORTAS-RS*DS*FACS  :  GOSUB  5070  : 

AS  C I , J 1 3 -SORTAS : CT-CT* 1 : R1 S  C  CT  3 -SORTAS : R2  C  CT  3 - 1 : R3C  CT  3 - J 1 
4000  ASCI , Jl*13-"000000" 

4010  ASC I 1 , Jl*l 3-"000000” 

4020  OPERS C FAC 3 -"CUR" 

4030  PL I NSC FAC 3 -RS 

4040  FOR  12-11*1  TO  1-1 

4050  ASCI2.J13-"-” 

4060  FOR  J2-J+1  TO  Jl-1 

4070  ASC 12 , J23-FACS 


4080  NEXT  J2 

4080  NEXT  12 
4100  RETURN 

4110  ’>>>  <<<  BOX  LEFT-DOWN 

4120  IF  DFLAG2-1  THEN  Il-D-1  ELSE  Il-I+C INTCABSC I-D3/233 

4130  IF  LFLAG2-1  THEN  Jl-L+1  ELSE  Jl-J-C INTCABSC J-L3/233 

4140  LS-ASCI+1, Jl-13 

41S0  US-ASCI-1, Jl  +  13 

4160  RS-ASCI+1, J+13 

4170  SORTAS- ”AA "+LS+F ACS  :  GOSUB  5070  :  ASC I 1 , J1 3-SORTAS 

4180  SORTAS-RS+LS+FACS  :  GOSUB  S070  : 

A* C  1 1 ,  J 1  -SORTAS :  CT-CT+1 :  R1SCCT 3-SORTAS:  R2CCT3-I1 R3CCT3-J 
4190  SORTAS-US+RS+FACS  :  GOSUB  5070  : 

AS  C I , J 1 -SORTAS : CT-CT+1 : R1SCCT3-S0RTAS: R2CCT3-I : R3CCT3-J 
4200  SORTAS-US+LS+FACS  :  GOSUB  5070  : 

AS  C I , J 1 3 -SORTAS : CT-CT+1 : R1SCCT3-S0RTAS: R2CCT3-I : R3CCT3-J1 
4210  ASC I 1+1, J3-”000000" 

4220  ASCI , Jl-13-"000000” 

4230  0PERSCFAC3-"BLD" 

4240  PL I NSC FAC 3 -LS 

4250  FOR  I 2-1+1  TO  11-1 
4260  ASC 12, J13-“~” 

4270  NEXT  12 

4280  FOR  J2-J1+1  TO  J-l 

4290  ASCII, J23-"-” 

4300  FOR  I 2-1+1  TO  11-1 

4310  ASC 12, J23-FACS 

4320  NEXT  12 

4330  NEXT  J2 

4340  RETURN 

4350  ■»>  <<<  BOX  LEFT-UP 

4360  IF  UFLAG2-1  THEN  Il-U+1  ELSE  Il-I-C INTCABSC I-U3/23 3 

4370  IF  LFLAG2-1  THEN  Jl-L+1  ELSE  Jl-J-C INTCABSC J-L 3/2 3  3 

4380  LS-ASC I-1.J1-13 
4390  DS-ASCI+1, Jl+13 

4400  RS-ASCI-1, J+13 

4410  SORTAS- ”DD"+LS+FACS  :  GOSUB  5070  :  ASC I 1 , J1 3-SORTAS 

4420  SORTAS-RS+LS+FACS  :  GOSUB  5070  : 

AS  C 1 1 , J  3-SORTAS : CT-CT+1 : R1 S  C  CT  3 -SORTAS : R2  C  CT 3- 1 1 : R3 C  CT  3 - J 
4430  SORTAS-DS+RS+FACS  :  GOSUB  5070  : 

AS  C I , J  3 -SORTAS :CT-CT+1:R1SCCT  3-SORTAS : R2  C  CT  3 - 1 : R3C  CT  3 - J 
4440  SORTAS-DS+LS+FACS  :  GOSUB  5070  : 

ASC I , J1 3-SORTAS: CT-CT+1 : R1S CCT 3-SORTAS : R2CCT 3- I : R3CCT3-J1 
4450  AS  CI1-1.J3- "000000 " 

4460  ASCI , J1 -13 -“000000” 

4470  OPERS  C  FAC  3- "BLU ” 

4480  PL I NSC FAC 3 -LS 

4490  FOR  I 2-1 1+1  TO  1-1 
4500  ASCI2.J13-"-" 

4510  NEXT  12 

4520  FOR  JZ-J1+1  TO  J-l 
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4530  ASCII, J23-"-" 

4540  FOR  12-11*1  TO  1-1 

4550  ASCI 2. J23-FACS 

4560  NEXT  12 

4570  NEXT  J2 

4580  RETURN 

4590  ’>>>  <<<  BOX  RIGHT-DOWN 

4600  IF  DFLAG2-1  THEN  Il-D-1  ELSE  I 1-I  +  C INTC ABSC I-D3/23 3 

4610  IF  RFLAG2-1  THEN  Jl-R-1  ELSE  Jl-J+C 1NTCABSC J-R3/233 

4620  LS-ASC1+1.J-13 

4630  US-ASCI-1, Jl-13 

4640  RS-ASCI+1, Jl*13 

4650  SORTAS-"BB“+RS*FACS  :  GOSUB  5070  :  ASC I 1 , J1 3-SORTAS 
4660  SORTAS-RS*LS+FACS  :  GOSUB  5070  : 

ASCII, J3-S0RTAS:CT-CT*1:R1SCCT3-S0RTAS:R2CCT3-I1:R3CCT3-J 
4670  SORTAS-US+LS+FACS  :  GOSUB  5070  : 

ASC I , J3-S0RTAS: CT-CT+1 : R1SCCT 3-SORT AS: R2CCT3-I : R3CCT3-J 
4680  SORTAS-US+RS+FACS  :  GOSUB  5070  : 

AS  C I , J 1 3 -SORTAS : CT-CT+ 1 : R1 SC  CT  3 “SDRTAS : R2  C  CT  3 - 1 : R3  C  CT  3 - J 1 
4690  ASCI1+1, J3-”000000” 

4700  AS  Cl ,  J 1  ♦  1  3  - ”000000 ” 

4710  0PERSCFAC3--BR0” 

4720  PLINSCFAC3-RS 

4730  FOR  I2-I+1  TO  I 1-1 

4740  ASC 12, J1 3““-“ 

4750  NEXT  12 

4760  FOR  J2-J*l  TO  Jl-1 

4770  ASCII, J23-”-” 

47B0  FOR  I 2-1*1  TO  I 1-1 

4790  ASC I 2, J2 3 -FACS 

4800  NEXT  12 

4810  NEXT  J2 

4820  RETURN 

4830  ■>>>  <<<  BOX  RIGHT-UP 

4840  IF  UFLAG2-1  THEN  Il-U+1  ELSE  I 1-I-C INTC ABSC I-U3 /23 3 

4850  IF  RFLAG2* 1  THEN  Jl-R-1  ELSE  J1-J*C INTCABSC J-R3/23 3 

4860  LS-ASC I -1 , J-l 3 

4870  DS-ASC 1*1, Jl-13 

4880  RS-ASC I-1.J1+13 

4890  SORTAS- ”CC"*RS*F ACS  :  GOSUB  5070  :  ASC  1 1 ,  J1  3-SORTAS 

4900  SQRTAS-RS*LS*FACS  :  GOSUB  5070  : 

AS C 1 1 , J  3-SORTAS : CT-CT+1 : R IS C  CT  3 -SORTAS : R2 C CT  3 - 1 1 : R3C CT  3 - J 
4910  SQRTAS-DS+LS+FACS  :  GOSUB  5070  : 

ASC I , J 3-SORTAS :CT-CT*1 :R1S CCT 3-SORTAS :R2CCT 3-1 :R3CCT3-J 
4920  SORTAS-DS*RS*FACS  :  GOSUB  5070  : 

ASC I , J1 3-SORTAS: CT-CT+1 : R1SCCT3-SDRTAS: R2CCT3-I : R3CCT3-J1 
4930  ASC 1 1-1, J 3- ”000000" 

4940  ASC I , Jl* 13- "000000" 

4950  0PERSCFAC3-"BRU” 

4960  PLINSCFAC3-RS 

4970  FOR  12-11*1  TO  1-1 
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4980 

4930 

5000 

5010 

5020 

5030 

5040 

5050 

5050 

5070 

5080 

5090 

5100 

5110 

5120 

5130 

5140 

5150 

5160 

5170 

5180 

5190 

5200 

5210 

5220 

5230 

5240 

5250 

5250 

5270 

5280 

5290 

5300 

5310 

5320 

5330 

5340 

5350 

5360 

5370 

5380 

5390 

5400 

5410 

5420 

5430 

5440 


ASCI2, J13-"-" 

NEXT  12 

FOR  J2-J+1  TO  Jl-1 

ASC I 1 , J23-”-” 

FOR  I 2-11+1  TO  1-1 

ASCI2, J23-FACS 

NEXT  12 

NEXT  J2 
RETURN 

■>>>  <<<  SORT  ROUTINE 

’SORTAS  IS  THE  ELEMENT  TO  BE  PUT  IN  NUMERICAL  ORDER 
N1S-M I DSC SORTAS ,1,23 
N2S-M I DS  C  SORTAS ,3,23 
N3S-M I D$  C  SORT AS , 5 , 2  3 

IF  N1S-”AA"  OR  N1S-’’BB"  OR  N1S-"CC”  OR  N1S-“DD"  GOTO  5270 

NC13-UALCN1S3 

NC23-UALCN2S3 

NC33-UALCN3S3 

FOR  X-l  TO  2 

FOR  Y-X+l  TO  3 

IF  NCX3  <-  NCY3  GOTO  5220 

H-NCX3 

NCX3-NCY3 

NCY3-H 

NEXT  Y 

NEXT  X 

IF  NC 1 3< 10  THEN  NlS-”0“+RIGHTSCSTRSCNC 1 3 3 , 1 3  ELSE 

N1S-RIGHTSCSTRSCNC 133,23 

GOTO  5330 

GOTO  5360 

NC23-UALCN2S3 

NC33-UALCN3S3 

IF  NC23  <-  NC33  GOTO  5330 

H-NC23 

NC23-NC33 

NC33-H 

IF  NC23< 10  THEN  N2S-”0"+RIGHTSCSTRSCNC23 3 , 1 3  ELSE 
N2S-RIGHTSCSTRSCNC233 ,23 

IF  NC  33< 10  THEN  N3S-"0"+R IGHTSCSTRSCNC33 3  ,  1 3  ELSE 

N3S-RIGHTSCSTRSCNC333 ,23 

50RTAS-N1S+N2S+N3S 

RETURN 

•>>>  <<<  PRINT  INSERTION  ORDER 

LPRINT  ORDERC23; AREAC23 ; 0RDERC33 ; "s " ; AREAC 33 ; ORDERC43 ; " : 
FOR  1-5  TO  FAC-1 

LPRINT  ORDERC I  3 ; AREAC 13;”  ” ; OPERSC I  3 ; "  ” ; TR I ANGLES C I  3 ; 
PL INC  I  3 

NEXT  I 
RETURN 

•>>>  <<<  CLEAR  FIRST  LINE 
LOCATE  1,1 


AREAC43 


SB 


5150 

PRINT  BLS 

5160 

LOCATE  1,1 

5170 

RETURN 

5180 

■>>>  <<<  AREA  CALCULATIONS 

5130 

DIN  ULI C503  , ULJCS03 , UR It 501 , URJC503 , LLI C50D , LLJC503 ,LRI C503 , LRJC503 

5S00 

CLS 

5510 

FAC-FAC-1 

sseo 

PLINC33-2 

5530 

PLINC13-3 

5510 

FOR  1-FAC  TO  2  STEP  -1 

5550 

AREAINCPLINC I 3  3-AREAINCPLINC I 3  3+AREAINC I 3 

5560 

NEXT  I 

S570 

AREATDT-AREAINC23 

S5B0 

ULI C  23-0: ULJ  C23-0 

5530 

URIC23-0:URJC23-1 

5600 

LLI C23-1 : LLJC23-0 

5610 

LRIC23-1: LRJC23-1 

5630 

DRW-2 

5630 

GOSUB  8550 

5610 

CARUE-3  : 1-2 

5650 

GOSUB  6110 

5660 

CARUE-1  : 1-3 

5670 

GOSUB  6530 

S6B0 

FOR  1-3  TO  FAC 

5630 

CARUE-0 

5700 

BOX 1-0 

5710 

B0X2-0 

5720 

FOR  11- I  TO  FAC 

5730 

IF  PLINC I 1 3-1  AND  LEFTSCOPERSC I 1 3 , 1 3-"C” 

THEN  CARUE' 

5710 

IF  PLINC 113-1  AND  LEFTSCOPERSC I 1 3 , 1 3-"B” 
B0X1-I1 : GOTO  5770 

THEN 

S7S0 

NEXT  11 

5760 

GOTO  5810 

5770 

FOR  I2-I1+1  TO  FAC 

57B0 

IF  PLINC I23-I  AND  LEFTSCDPERSC 123 , 13-"C” 

THEN  CARUE' 

5730 

IF  PLINC  I23-I  AND  LEFTSCOPERSC  123 , 1 3-”B’’ 
B0X2-I2: GOTO  5810 

THEN 

5B00 

NEXT  12 

5B10 

•CONTINUE 

5820 

IF  AREACB0X23>AREACB0X1 3  THEN  SWAP  B0X1.B0X2 

5830 

IF  CARUE-0  GOTO  5880 

5810 

IF  OPERS C CARUE 3 -"CLU"  OR  0PERSCCARUE3-”CRU”  THEN 
6110  : GOTO  5BB0 

GOSUB 

S850 

IF  OPERS C CARUE 3 - ”CLD ”  OR  OPERS C CARUE 3- ”CRD ”  THEN 
: GOTO  5880 

GOSUB  6270 

5860 

IF  OPERS C CARUE 3- "CDR "  OR  OPERS C CARUE 3 -”CUR”  THEN 
: GOTO  5880 

GOSUB  6130 

5870 

IF  OPERS C CARUE 3 -"CDL”  OR  OPERS C CARUE 3 -"CUL"  THEN 

GOSUB  6530 

5880 

IF  BOXl-O  GOTO  6060 

5830 

B0X-B0X1 

5300 

IF  0PERSCB0X1 3-"BLD"  THEN  GOSUB  6750: GOTO  5310 

S3 


5910 

IF 

OPERSC  B0X1 3- “BLU” 

THEN 

GOSUB 

6930 : GOTO 

5940 

5930 

IF 

0PERSCB0X1 3-”BRD” 

THEN 

GOSUB 

7090: GOTO 

5940 

5930 

IF 

0PERSCB0X1 3-“BRU” 

THEN 

GOSUB 

7360 

5940 

IF 

B0X3-0  GOTO  6060 

5950 

FLAGFIT-0 

5960 

IF 

OPERS  C  B0X2  3  -  ”BLD  ’’ 

THEN 

GOSUB 

7430: GOTO 

6000 

5970 

IF 

OPERSC B0X3 3 -“BLU” 

THEN 

GOSUB 

7710: GOTO 

6000 

5980 

IF 

OPERSCBOX33-”BRD" 

THEN 

GOSUB 

7990: GOTO 

6000 

5990 

IF 

OPERSCBOX33-“BRU" 

THEN 

GOSUB 

8370 

6000 

IF 

FLAGFIT-1  GOTD  6060 

6010 

B0X-B0X3 

6030 

IF 

OPERS  C  B0X3 3 “ "BLD ” 

THEN 

GOSUB 

6750: GOTO 

6060 

6030 

IF 

OPERS  C  B0X3 3 - "BLU ” 

THEN 

GOSUB 

6930: GOTO 

6060 

6040 

IF 

OPERS  C  B0X2  3 - "BRD” 

THEN 

GOSUB 

7090: GOTO 

6060 

6050 

IF 

OPERSC BOX33-“BRU” 

THEN 

GOSUB 

7360 

6060  NEXT  I 

6070  LOCATE  S3 , X 

6080  INPUT  "Uould  you  like  a  list  oF  coordinates  printed  CY/N3" ;ANSS 

6090  IF  ANSS--Y"  OR  ANSS-"y”  THEN  GOSUB  8580 
6100  RETURN 

6110  '>>>  <<<  CLU  OR  CRU 

6130  0 1 STUP-AREA I N  C  CARUE 3 /AREA INC  I 3*CLLI C 1 3-ULI t I  3  3 

6130  LLICCARUE3-LLICI3 

6140  LRICCARUE3-LRICI3 

6150  LLICIl-LLIcn-DISTUP 

6160  LRI C I 3-LR1 C 1 3-DISTUP 

6170  ULI CCARUE3-LLI C I  3 

6180  URICCARUE3-LRIC13 

6190  ULJCCARUE3-ULJCI3 

6300  UR J  C  CARUE  3 -UR J  C I  3 

6310  LLJCCARUE3-LLJCI3 

6330  LR J  C  CARUE  3 -LR J 1 1  3 

6330  DRUI-CARUE 

6340  GOSUB  8550 

6350  AREA  INCH  -AREA  I N  C I  3 - AREA I NC  CARUE  3 

6360  RETURN 

6370  ’>>>  <<<  CLD  OR  CRD 

6380  D I STDUIN-AREA I N  C  CARUE  3 /AREA I N  C I  3 • C  LL I C I  3 -UL I C I  3  3 

6390  UL I C  CARUE  3 -UL I C I  3 

6300  URI CCARUE3-URI C I  3 

6310  UL  I C  I  3  -UL  l C I  3 +D I STDUIN 

6330  URICI3-URICI3+D1STDUIN 

6330  LL 1 C  CARUE  3 -UL 1 C 1 3 

6340  LRI CCARUE3-URI C I  3 

6350  UL J  t  CARUE  3 -UL J  C I  3 

6360  UR J  C  CARUE  3 -UR J  C I  3 

6370  LL J  C  CARUE  3 -LL J 1 1  3 

6380  LR J  C  CARUE  3 -LR J  C I  3 

6390  DRUI-CARUE 

6400  GOSUB  B550 

6410  AREAINC I 3-AREAINC I  3 -AREA INC  CARUE  3 


6420  RETURN 

6*130  *>>>  <<<  COR  OR  CUR 

6440  DISTRT-AREAINCCARUE 3 /AREAINC I 3-CURJC I D-ULJ  C  I  3  D 

6450  ULJCCARUE3-ULJCI3 

6460  LL  J  C  CARUE  3 -LL J  C I  3 

6470  ULJCI3-ULJCI3+OISTRT 

6480  LLJC  I  3-LLJC  I  D+'OISTRT 

6430  UR J  C  CARUE  3 -UL J  C I  3 

6500  LR J  C  CARUE  3 -LL J  C I  3 

6510  ULI CCARUE3-ULI C I  3 

6520  URIC  CARUE  3 -UR I C I  3 

6530  LLI CCARUE3-LLI C I  3 

6540  LR I C  CARUE  3 -LR I C 1 3 

6550  ORU-CARUE 

6560  GOSUB  8550 

6570  AREA INC  I  3 “AREA INC  I  3 -AREA INC CARUE 3 

6580  RETURN 

6590  ’>>>  <<<  COL  OR  CUL 

6600  0ISTLT-AREAINCCARUE3/AREAINC I 3»CURJC I 3-ULJC I  3  3 

6610  URJCCARUE3-URJC I  3 

6620  LR J  C  CARUE  3 -LR J  C 1 3 

6630  URJC I 3-URJC I 3-OISTLT 

6640  LRJC I 3-LRJC I 3-OISTLT 

6650  UL J  C  CARUE  3 -UR J  C I  3 

6660  LLJCCARUE3-LRJC I 3 

6670  UL I C  CARUE  3 -UL I C I  3 

6680  URI CCARUE3-URI C I  3 

6690  LLI C CARUE 3 -LLI C I  3 

6700  LR I C  CARUE  3 -LR I C I  3 

6710  DRU-CARUE 

6720  GOSUB  8550 

6730  AREAINC I 3-AREAINC I  3 -AREA INC CARUE 3 

6740  RETURN 

6750  ’>>>  <<<  BLD 

6760  0 1 STLT- C  AREA INC  BOX  3 /AREA I N C I  3  3  A . 5*  C  UR J  C I  3 -UL J  C 1 3  3 

6770  DISTDWN-CAREAINCB0X3/AREAINC I33A.5*CLRICI 3-URI C 133 

6780  URI CB0X3-UFI C I  3 

6790  URJCBOX3-URJCI3 

6800  URI C I 3-URI C I 3+D15TDWN 

6810  URJC I 3-URJC I 3-DISTLT 

6820  LLI CB0X3-URI C 1 3 

6830  LLJCBOX3-URJCI3 

6840  ULI CB0X3-ULI C 1 3 

6850  UL J  C  BOX  3 -UR J  C I  3 

6860  LRI CB0X3-URI C I  3 

6870  LRJCB0X3-LRJC I  3 

6880  DRW- BOX 

6890  GOSUB  8550 

6300  AREAINC 1 3 -AREAINC I  3 -AREA INC BOX 3 

6910  RETURN 

6920  ’>>>  <<<  BLU 
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6330  DISTLT-CAREAINCB0X3 /AREAINC I  3 3A .S»CLRJC I 3-LLJCI 3  3 

6940  DISTUP-C AREAINCB0X3/AREAINC I  3  3  A  . 5* C  LR I C 1 3-URI C I  3  3 

6950  LRICBOX3-LRIC13 

6960  LRJCBOX3-LRJCI3 

6970  LRICI3-LRICI3-DISTUP 

6980  LR J  C I  3 -LR J 1 1  3 -0 1 STLT 

6990  ULI CB0X3-LRI C 1 3 

7000  UL J  C  BOX  3 -LR J  C I  3 

7010  LLI CB0X3-LL1 C 1 3 

7020  LLJCB0X3-LRJC I  3 

7030  URI CB0X3-LRI C I  3 

7040  UR J  t  BOX  3 -UR J  C I  3 

7050  ORU-BOX 

7060  GOSUB  9550 

7070  AREA1NC I 3-AREAINC I 3-AREAINCBOX3 

7000  RETURN 

7090  ’>>>  <<<  BRD 

7100  DISTRT-C AREAINC BOX 3 /AREAINC I  3  3A .5»CURJC 1 3-ULJC 1 3  3 

7110  D I STDUN-C AREA INC  BOX  3 /AREAINC I33A.5*CLLICI 3-UL I C I  3  3 

7180  ULI CBQX3-UL! C  I  3 

7130  ULJCBOX3-ULJCI3 

7140  ULI C I 3-ULI C I 3+DISTDUN 

7150  ULJC I 3-ULJC I 3+DISTRT 

7160  LRI [BOX 3-ULI C I  3 

7170  LRJC BOX 3-ULJC 13 

71 BO  URI CB0X3-URI C I  3 

7190  UR J  C  BOX 1 -UL J  C I  3 

7200  LL I C  BOX  3 -UL I C I  3 

7210  LL J  C  BOX  3 -LL J  C I  3 

7220  DRU-BOX 

7230  GOSUB  8550 

7240  AREAINC I 3-AREAINC I  3 -AREA INC BOX 3 

7250  RETURN 

7260  *>>>  <<<  BRU 

7270  0 1 STRT- C  AREA INC  BOX  3 /AREA I NC 1 3  3  A . 5* C  LR J  C I  3 -LL J C  I  3  3 

7280  0 1 STUP- C AREA INC B0X3/ AREA INC  I  3  3 A . 5»CLLI C I 3-ULI C 133 

7290  LLI CB0X3-LLI C I  3 

7300  LLJC BOX 3-LLJCI 3 

7310  LLI C I 3-LL1 C I 3-OISTUP 

7320  LLJC I 3-LLJC I 3+DISTRT 

7330  URI CB0X3-LLI C I  3 

7340  UR J  C  BOX  3 -LL J  C 1 3 

7350  LRICBOX3-LRICI3 

7360  LRJCB0X3-LLJC 1 3 

7370  ULICBOX3-LL1CI3 

7380  UL J  C  BOX  3 -UL J  C I  3 

7390  ORU-BOX 

7400  GOSUB  B550 

7410  AREAINC I 3-AREAINC I  3-AREA INC BOX 3 

7420  RETURN 

7430  ’>>>  <<<  BLD  CORRECTIONS 


7440  IF  0PERSCB0X13-”BLU”  GOTO  7580  ’IF  NOT  IT”S  BRO 
7450  LLJCI3-ULJCI3 

7460  uLicn-uRicn 

7470  AREAINC 1 3-CLRJC I 3-LLJC 1 3 3*CLLI C I 3-ULI C I  3  3*AREAT0T 

7480  IF  AREAINCI3*.9B<AREAINCBOX23  THEN  DRU-BQX1 : GOSUB  BB60:G0SUB  8700  ELSE 
GOTO  7700 

7430  LLJC I 3-ULJC I  3 
7500  URJCBOX13-ULJCI3 

7510  LLI CB0X1 3-CAREAINCB0X1 3 /AREATDT/CURJC B0X1 3-ULJC B0X1 3  3  3+ULI CB0X1 3 

7520  LRJCBOX13-URJCBOX13 

7S30  LRI CB0X13-LLI CB0X13 

7540  0RU-B0X1 

7550  GOSUB  0550 

7580  LINE  C I  NT  C  UR  JC  BOX  1 3  *400  3 -•'GO,  INTCURI C BOX 1 3*1503+253  - 
C INTCURJC I  3 *400 3 ♦60, INTtURI 1 1  3*1503*253 
7570  GOTO  7700 

7580  LRJCI3-URJCI3 
7590  LLICI3-LRICI3 

7600  AREAINC  I  3-CLRJC  1 3 -LLJC  1 3  3*CLLI  C  I  3-ULI  C 1 3  3*AREATDT 

7610  IF  AREAINC I  3* .9B< AREAINCBOX23  THEN  DRW-B0X1 : GOSUB  BB60:G05UB  8740  ELSE 

GOTO  7700 

7620  LRICI3-LLICI3 
7630  URICBOX13-LLICI3 

7640  LL J  C  BOX 1 3  —  C  ARE A I N  C  BOX 1 3 / AREATOT / C  LR I C  BOX 1 3 -UR I C  BOX 1 3  3  3  *LR J  C  BOX 1 3 

7650  ULICBOXll-URKBOXn 

7660  ULJCBOX13-LLJCBOX13 

7670  DRW-B0X1 

7680  GOSUB  8550 

7690  LINE  C I NT C UR JC BOX 1 3*4003*60, INTCURI CB0X1 3*1503*253  - 
C INTCURJC I  3*4003*60, INTCURI C I  3*1503*253 
7700  RETURN 

7710  ’>>>  <<<  BLU  CORRECTIONS 

7720  IF  0PERSCBQX1 3“"BRU"  GOTO  7B60  ’IF  NOT  IT”5  BLD 
7730  URJC 1 3-LRJC 1 3 

7740  ULICI3-URICI3 

7750  AREA I NC I  3 - C  LR J C I  3 -LL J C I  3  3*  C LL I C 1 3 -UL I C 1 3  3 *AREATOT 

7760  IF  AREAINC I  3* ,9B< AREAINC BOX23  THEN  0RW-B0X1 : GOSUB  8860: GOSUB  BB20  ELSE 
GOTO  7980 

7770  URICI3-ULICI3 
7700  LRICBOX13-ULICI3 

7790  ULJCB0X13— C AREAINCB0X1 3/AREATOT/CLRI CB0X1 3-URI CB0X1 3 3 3*URJCB0X1 3 

7800  LLI CB0X1 3-LRI CB0X1 3 

7810  LL J  C  BOX 1 3 -UL  J C  BOX 1 3 

7820  DRW” BOX 1 

7830  GOSUB  8550 

7840  LINE  C 1NTCLRJCBQX1 3*4003*60 , INTCLRI CB0X1 3*1503+253  - 
C INTCLRJC I  3*4003+60, INTCLRI C I  3*1503+253 
7850  GOTO  79B0 
7860  ULJCI3-LLJCI3 
7070  LLICI3-LRICI3 

7880  AREAINC I 3-CLRJC I 3-LLJC I  3  3*CLL I C I 3-ULI C I  3  3 "AREATOT 
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7830  IF  AREA INC  I la . 98< AREA  I N[ 80X21  THEN  DRU-B0X1 : GOSUB  BB60: GOSUB  8700  ELSE 
GOTO  7980 

7300  LLJCI1-ULJCI1 

7310  LRJCB0X11-ULJCI1 

7320  ULI CB0X1 1--C AREA INC B0X1 1/ AREATOT/ C LRJC B0X1 1-LLJC  BOX 1 1 1 1+LL I CB0X1 1 

7930  URJCB0X1 1-LRJCB0X1 1 

79*10  URICB0X1 1-ULI CB0X11 

7950  DRW- BOX 1 

7960  GOSUB  8550 

7970  LINE  C INTCLRJC B0X1 3 *400 3 *60 , INTCLRI CB0X1 1*1501+251  - 
C INTCLRJC 1 1*4001+60 , INTCLRI C 1 1*1501+251 
7980  RETURN 

7990  ’>>>  <<<  BRO  CORRECTIONS 

8000  IF  0PERSCB0X11-"BL0"  GOTO  8140  ’IF  NOT  IT”S  BRU 

8010  LLJCI1-ULJCI1 

8020  LRICI1-LLICI1 

8030  AREAINC 1 1-CLRJC 1 1-LLJC 1 1 1*CLLI C 1 1-ULI C 1 1 1*AREAT0T 

8040  IF  AREAINC 1 1* . 9B< AREAI NC B0X21  THEN  DRU-B0X1 : GOSUB  8860: GOSUB  8740  ELSE 
GOTO  8260 

8050  LRICI1-LLICI1 

8060  ULICB0X11-LLICI1 

8070  LR J  C  BOX 1 1 - C  AREA I N  C  BOX 1 1 / AREATOT / C  LL I C  BOX 1 1 -UL I C  BOX 1 1 1 1 +LL J  C  BOX 1 1 

8080  URICB0X1 1-ULI CB0X11 

8090  URJCB0X1 1-LRJCB0X1 1 

8100  0RU-BQX1 

8110  GOSUB  8550 

8120  LINE  C I NT CULJC BOX 1 1*4001+60, I NT C ULI CB0X1 1*1501+251  - 
C I NT CULJC 1 1*4001+60, INTCULI C 1 1*1501+251 
8130  GOTO  8260 

8140  LRJC 1 1-URJC 1 1 

8150  URI C I 1-ULI C I 1 

8160  AREAINC 1 1-CLRJC 1 1-LLJC 1 1 1*CLLI C 1 1-ULI C 1 1 1*AREAT0T 

8170  IF  AREAINC I 1*.38<AREAINCB0X21  THEN  0RW-B0X1 : GOSUB  BB60:G0SUB  B7B0  ELSE 
GOTO  B260 

8180  LRJC I 1-URJC II 

8190  ULJCB0X1 1-URJC II 

8200  LRI CB0X1 1-C AREA INC B0X1 1 /AREATOT/ CUR JCB0X1 1-ULJCB0X1 1 1 1+URI CB0X1 1 

8210  LLJCB0X1 1-ULJC  B0X1 1 

B220  LLI CB0X1 1-LRI CB0X1 1 

8230  DRW-B0X1 

8240  GOSUB  BS50 

8250  LINE  C INTCULJC B0X1 1*4001+60 , INTCULI C BOX 1 1*1501+251  - 
C INTCULJC 1 1*4001+60, INTCULI C 1 1*1501+251 
8260  RETURN 

8270  ’>>>  <<<  BRU  CORRECTIONS 

8280  IF  OPERS C BOX 1 1 “ "BRO "  GOTO  8420  ’IF  NOT  IT”S  BLU 
8290  URJC 1 1-LRJC 1 1 

8300  LRI C 1 1-LLI C 1 1 

8310  AREAINC 1 1-CLRJC I 1-LLJC I 11*CLLI C I 1-ULI C 1 1 1 -AREATOT 

8320  IF  AREAINC I 1*.98<AREAINCB0X21  THEN  DRW-B0X1 : GOSUB  8860: GOSUB  8780  ELSE 

GOTO  8540 


i 
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6330  LRJCI3-URJCI3 

0340  LLJCBOXn-URJCn 

8350  UR I C  B0X1 3  —  C AREA I N t  BOX 1 3 / AREATOT/ C  LR JC BOX 1 3 -LL J C  BOX 1  3  3  3 +LR 1 C  BOX 1 3 

8360  UL J  C  BOX 1 3 -LL J  C  BOX 1 3 

8370  ULICBOX13-URICBOX13 

8380  DRW-B0X1 

8390  GOSUB  B550 

8400  LINE  C I NTCLLJCBQX1 3*4003+60 , INTCLLI C 80X13 *150 3+25 3  - 
C INTCLLJC I  3 •4003+60, INTCLLI C I 3*1503+253 
8410  GOTO  8540 

8420  ULJCI3-LLJCI3 

8430  URItn-ULUI3 

8440  AREAINC I 3-CLRJC I 3-LLJt 13  3*CLLI C I 3-ULI C I  3  3 “AREATOT 

8450  IF  AREAINC 13* ,98< AREAINC  B0X2 3  THEN  DRW-B0X1 : GOSUB  8B60:G0SUB  8820  ELSE 
GOTO  8540 

8460  URICI3-ULICI3 

8470  LLICB0X1 3-ULI C 13 

B4B0  URJCB0X1 3-CAREAINCB0X1 3/AREATQT/CLLI CB0X1 3-ULI CB0X13  3  3+ULJCB0X1 3 

8490  LR1 CB0X1 3-LLI CB0X1 3 

8500  LRJCB0X1 3-URJCB0X1 3 

8510  0RU-B0X1 

8520  GOSUB  8SS0 

0530  LINE  C I NT CLLJCB0X1 3*4003+60, INTCLLI CBQX1 3*1503+2S3  - 
C INTCLLJC I  3*4003+60, INTCLLI C I  3*1503+253 
0540  RETURN 
0550  •>>>  <<<  BOX 

8560  L I NE  C I NT  C  UL J  C  DRW  3  *400  3  +60 , I NT  C  UL I C  DRW  3*1 SO  3  +25  3  - 

C I NT  C  LR J  C  DRW  3  *400  3  *60 , INTCLRI CDRW3*1503+253 , ,B 
8570  RETURN 

0500  •>>>  «<  PRINT  COORDINATES 


0590 

FOR  1-2  TO  FAC 

8600 

LPRINT 

ORDERC I 3 

8610 

LPRINT 

USING 

”#»# . «««« ” j  UL I C I 3 ; UL J C I 3 

8620 

LPRINT 

m 

H  . 

• 

8630 

LPRINT 

USING 

UR I C I 3 ; URJC I 3 

8640 

LPRINT 

USING 

; LLI C I 3 ; LLJC I 3 

B650 

LPRINT 

N 

H  . 

• 

8660 

LPRINT 

USING 

”#*# . «««« " ( LR I C I 3 ; LRJ C I 3 

8670 

LPRINT 

8680 

NEXT  I 

8690 

RETURN 

8700 

•>>>  <<<  PUSH  LEFT 

8710 

ULJC I 3-ULJC I 3- 

.01 

8720  AREAINC I 3-CURJC I 3-ULJC I  3  3*CLLI C I 3-ULI C I  3  3*AREAT0T 
8730  IF  AREAINC I  3* . 9B< AREAINCBOX23  THEN  GOTO  8710: ELSE  RETURN 
8740  ’>>>  <<<  PUSH  DOWN 

8750  LLI C I 3-LLI C I  3+ .01 

0760  AREAINC I 3-CURJC I 3-ULJC I  3  3*CLLI C I 3-ULI C 1 3  3 "AREATOT 
8770  IF  AREAINC I  3* ,98<AREAINCB0X23  THEN  GOTO  8750: ELSE  RETURN 
8780  ’>>>  <<<  PUSH  RIGHT 

8790  URJC I 3-URJC I  3+ .01 


8800  AREA  INC  I  3“CURJC  I  3-ULJC  I  3  3  *  C  l_l_  I C  I  3-ULI C  1 3  3*AREAT0T 
8810  IF  AREAINC 1 1* .98<AREAINCB0X23  THEN  60T0  8790: ELSE  RETURN 
8820  ’>>>  <<<  PUSH  UP 

8830  ULI Cl 3-ULI Cl  3-. 01 

8840  AREAINC I 3-CURJC I 3-ULJC 1 3  3*CLLI C I 3-ULI C I  3  3*AREAT0T 
8850  IF  AREAINC I  3* .9B< AREAINCBOX23  THEN  SOTO  8830: ELSE  RETURN 
8860  ’>>>  <<<  UNBOX 

8870  LINEC INTCULJCDRW3*4003+B0, INTCULICDRW3*1503+E53  - 
C 1 NT CLRJCDRU 3 *4003+60, INTCLRI CDRUI3*lS03+253 , 0. B 
8880  RETURN 
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APPENDIX  D 


OUTPUT  FROM  EXAMPLE  II 

Run 

Rendon  number  seed  C-3276B  to  3276']'''  2 

You  will  need  to  input  the  Filename  for  the  data  you  wont  to  use 
Would  you  like  o  list  of  files  on  the  disk  CY/NJ-*  N 
Enter  ony  filename  with  .DAT  for  on  extension7  3-12R 


If  you  need  an  X  value  other 
if  not  press  return.7 


than  -1024  enter  it  at  the  prompt. 


NUMBER 
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DELTAHEDRON  INSERTION  ORDER 


12  3 

4  5  6 

7 

B 

9  10  11 

12 

INSERTING 

UERTEX 

5 

IN 

TRIANGLE 

1 

2 

4  010204 

INSERTING 

UERTEX 

6 

IN 

triangle 

2 

3 

4  020304 

INSERTING 

UEPTEX 

7 

IN 

TRIANGLE 

5 

2 

4  050204 

INSERTING 

UERTEX 

e 

IN 

TRIANGLE 

2 

3 

6  020306 

INSERTING 

UERTEX 

9 

IN 

TRIANGLE 

2 

3 

B  02C30B 

INSERTING 

UERTEX 

10 

IN  TRIANGLE 

2 

3 

9  0203C9 

INSERTING 

UEPTEX 

u 

IN  TRIANGLE 

1 

2 

3  010203 

INSERTING 

UERTEX 

12 

IN  TRIANGLE 

1 

4 

5  010405 

TOTAL  DELTAHEDRON  ADJACENCY  SCORE  IS  614 
INCIDENCE  MATRIX: 


1 

2 

3 

4 

5 

6 
7 
B 
3 
10 
11 
12 


AUUU-----AA 
AA0EUEI  I  A  - 
a-a-aouA- 

E  E  U  - - U 

-  E  -  - A 


0  - - 

I  -  - 


E -ample  II  Dei to^adran  Heuristic  Output 


Vi  J 


:*1 


m 


1 


liiiiiiiiiiiiii 


11111111111 


2  - 
2  - 
8  - 
1  -  2  - 

-  e  - 
1  -  2  - 

2  - 
1  -  2  - 

2  0 
2  - 
2  - 
2  - 
1  -  2  - 

2  - 
1  -  2  - 

2  - 
2  - 
2  - 
2  - 


1  - 


555555555555555 


55S555S55555S55 
555555555555555 
555555555555  5*  55 
555555555555555 
555555555555555 
555555555555555 
555555555555555 
555555555555555 


77777777777777* 


7*T7777777777 


7  7  7  7  7 
"’  7  7  7  7 


7  7  7  7  7  - 


77777777777777 

777"’7777777777 

77777777777777 

77777777777777 


777777777 
77777  7  777 


-12121212121212121212121212 

-12121212121212121212121212 

-12121812121212121212121212 

-12121212121212121212121215 

-12121212121212121212121212 

-12121212121212121212121212 

-12121212121212121212121212 

-12121212121212121212121212 

-12121212121212121212121212 

-12121212121212121212121212 

-12121212121212121212121212 

-12121212121212121212121212 

-12121212121212121212121212 

-12121212121212121212121212 

-12121212121212121212121212 

-12121212121212121212121212 

-12121212121212121212121212 

-12121212121212121212121212 

-12121212121212121212121212 


1  - 
1  -  2 
1  -  2 
1  -  2 
1  -  2 
1  -  2 
1  -  2 
1  -  2 
1  -  2 
1  -  2 
1  -  2 
1  -  2 
1  -  2 
1  -  2 
1  -  2 
1  -  2 
1  -  2 
1  -  2 
1  -  2 
1  -  2 


-0B8BBBBBBB9BBB-3333 

-BBBBBBBSBBBBBB-3333 

-BBBBBBBBBBBB8B-3333 

0BBBBBB8BBBBBBB03333 

- 3333 

-99939939999399-3333 

09999999999999903333 

- 3333 

-1010101010101010101010101010  -  3  3  3  3 


3  3  3 
3  3  3 
3  3  3 
3  3  3 
3  3  3 
3  3  3 
3  3  3 
3  3  3 
3  3  3 


3  3  3  3  3  3 
3  3  3  3  3  3 


3  3  3  3  3  3 
3  3  3  3  3  3 


3  3  3  3  3  3 
3  3  3  3  3  3 


3  3  3  3  3  3 
3  3  3  3  3  3 


3  3  3  3  3 


3  3  3  3  3  3 
3  3  3  3  3  3 


333333333333333 


3  3  3  3  3  3 


1  -  2 
1  -  2 
1  -  2 
1  -  2 
1  -  2 
1  -  2 
1  -  2 
1  -  2 
1  -  2 


-111111111111111111111111111111111111111111111111111111111111 
-111111111111111111111111111111111111111111111111111111111111 
-111111111111111111111111111111111111111111111111111111111111 
-llllllllllllllllllllllllllllllllllllllllllllllllllllllllllll 
-111111111111111111111111111111111111111111111111111111111111 
-111111111111111111111111111111111111111111111111111111111111 
-111111111111111111111111111111111111111111111111111111111111 
-111111111111111111111111111111111111111111111111111111111111 
-111111111111111111111111111111111111111111111111111111111111 
-111111111111111111111111111111111111111111111111111111111111 
-111111111111111111111111111111111111111111111111111111111111 
-111111111111111111111111111111111111111111111111111111111111 
-111111111111111111111111111111111111111111111111111111111111 
-111111111111111111111111111111111111111111111111111111111111 
-111111111111111111111111111111111111111111111111111111111111 
-111111111111111111111111111111111111111111111111111111111111 
-111111111111111111111111111111111111111111111111111111111111 
-mu  liiiiiiiiiiiiiiiiiiiiiiiiiinini  miiiiiiiiiiiimii  li 


1111111111111111111111111111111111 


E  temple  II  Cc^densed  A*  nctnx 


Copy  available  to  DT1C  does  not 
ymH  hilly  legible  reproduction 


115 


B 


0.0000 

0.5640 

0 . 0000 

0 . 0000 

0.0000 

0.5848 

3 

0.5048 

1 .0000 

0.0563 

0.0563 

0.8976 

1 . 0000 

4 

0.5048 

1.0000 

0.5775 

0.5775 

0.5048 
1 .0000 

5 

0.7404 

0.8155 

0.0795 

0.6197 

0.5848 

0.8155 

6 

0 . 5840 
0.8976 

0.5014 

0.5014 

0.5048 

0.0976 

7 

0.5848 

0.7404 

0.6137 

0.6197 

0.5848 

0.7404 

8 

0.5848 

0.8975 

0 . 4466 
0.4466 

0.584B 

0.0976 

9 

0.5048 

0 . 8976 

0.3157 

0.3157 

0.5B48 
0 . 8976 

10 

0.5048 

0.8976 

0.1049 

0.1B49 

0.5848 
0 . 0976 

11 

0.5848 

1 .0000 

0.0000 

0 . 0000 

0.5848 
1 . 0000 

10 

0,0155 

1.0000 

0.6197 

0.6197 

0.0155 
1 . 0000 

1 .0000 
1 . 0000 


0 . 1643 
0.5775 


0.6197 

0.6197 


1 .0000 
1 .0000 


0.5775 

0.S77S 


0.0795 

0.8795 


0 . 5214 
0.5614 


0.4456 

0.4456 


0.3157 

0.3157 


0.0563 

0.0563 


1 .0000 
1 .0000 
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OUTPUT  FROM  EXAMPLE  III 

RUN 

Random  number  seed  C -33768  to  3376737  1 

You  will  need  to  input  the  filename  for  the  doto  you  wont  to  use. 
Would  you  like  a  list  of  files  on  the  disk  CY/N1?  n 
Enter  any  filename  with  .DAT  for  on  extension?  FOULDS 

If  you  need  an  X  value  other  than  -103*1  enter  it  ot  the  prompt, 
if  not  press  return.?  -1 
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INSERTING  UERTEX  7  IN  TRIANGLE  B  10  3  140718 

INSERTING  UERTEX  14  IN  TRIANGLE  80  19  18  11080B 

INSERTING  UERTEX  15  IN  TRIANGLE  IB  88  8  100305 

INSERTING  UERTEX  8  IN  TRIANGLE  8  10  B  050714 
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